Investigating Second Eigenvalue

  if (require("PageRank")) {
      library(PageRank)
    }else{
      devtools::install_github("ryangreenup/PageRank")
      library(PageRank)
    }

  library(pacman)
  pacman::p_load(PageRank, devtools, Matrix, igraph, mise, tidyverse, rgl, latex2exp)
#  mise()

Looking at Density

Constants

Define some constants

p    <- seq(from = 0.01, to = 0.99, length.out = 10)
beta <- seq(from = 1, to = 20, length.out = 4)
sz <- seq(from = 100, to = 10, length.out = 5)
input_var <- expand.grid("p" = p, "beta" = beta, "size" = sz)
input_var

Function to Build Graph

random_graph <- function(p, beta, size) {
      g1 <- igraph::erdos.renyi.game(n = size, p)
      A <- igraph::get.adjacency(g1) # Row to column
      A <- Matrix::t(A)

      A_dens <- mean(A)
      T      <- PageRank::power_walk_prob_trans(A, beta = beta)
      tr     <- sum(diag(T))
      e2     <- eigen(T, only.values = TRUE)$values[2] # R orders by descending magnitude
      return(c(abs(e2), mean(A), tr))
}

Return results

Map the function

nc <- length(random_graph(1, 1, 1))
Y <- matrix(ncol = nc, nrow = nrow(input_var))
for (i in 1:nrow(input_var)) {
  X <- as.vector(input_var[i,])
  Y[i,] <-  random_graph(X$p, X$beta, X$size)
  print(i/nrow(input_var))
}
[1] 0.005
[1] 0.01
[1] 0.015
[1] 0.02
[1] 0.025
[1] 0.03
[1] 0.035
[1] 0.04
[1] 0.045
[1] 0.05
[1] 0.055
[1] 0.06
[1] 0.065
[1] 0.07
[1] 0.075
[1] 0.08
[1] 0.085
[1] 0.09
[1] 0.095
[1] 0.1
[1] 0.105
[1] 0.11
[1] 0.115
[1] 0.12
[1] 0.125
[1] 0.13
[1] 0.135
[1] 0.14
[1] 0.145
[1] 0.15
[1] 0.155
[1] 0.16
[1] 0.165
[1] 0.17
[1] 0.175
[1] 0.18
[1] 0.185
[1] 0.19
[1] 0.195
[1] 0.2
[1] 0.205
[1] 0.21
[1] 0.215
[1] 0.22
[1] 0.225
[1] 0.23
[1] 0.235
[1] 0.24
[1] 0.245
[1] 0.25
[1] 0.255
[1] 0.26
[1] 0.265
[1] 0.27
[1] 0.275
[1] 0.28
[1] 0.285
[1] 0.29
[1] 0.295
[1] 0.3
[1] 0.305
[1] 0.31
[1] 0.315
[1] 0.32
[1] 0.325
[1] 0.33
[1] 0.335
[1] 0.34
[1] 0.345
[1] 0.35
[1] 0.355
[1] 0.36
[1] 0.365
[1] 0.37
[1] 0.375
[1] 0.38
[1] 0.385
[1] 0.39
[1] 0.395
[1] 0.4
[1] 0.405
[1] 0.41
[1] 0.415
[1] 0.42
[1] 0.425
[1] 0.43
[1] 0.435
[1] 0.44
[1] 0.445
[1] 0.45
[1] 0.455
[1] 0.46
[1] 0.465
[1] 0.47
[1] 0.475
[1] 0.48
[1] 0.485
[1] 0.49
[1] 0.495
[1] 0.5
[1] 0.505
[1] 0.51
[1] 0.515
[1] 0.52
[1] 0.525
[1] 0.53
[1] 0.535
[1] 0.54
[1] 0.545
[1] 0.55
[1] 0.555
[1] 0.56
[1] 0.565
[1] 0.57
[1] 0.575
[1] 0.58
[1] 0.585
[1] 0.59
[1] 0.595
[1] 0.6
[1] 0.605
[1] 0.61
[1] 0.615
[1] 0.62
[1] 0.625
[1] 0.63
[1] 0.635
[1] 0.64
[1] 0.645
[1] 0.65
[1] 0.655
[1] 0.66
[1] 0.665
[1] 0.67
[1] 0.675
[1] 0.68
[1] 0.685
[1] 0.69
[1] 0.695
[1] 0.7
[1] 0.705
[1] 0.71
[1] 0.715
[1] 0.72
[1] 0.725
[1] 0.73
[1] 0.735
[1] 0.74
[1] 0.745
[1] 0.75
[1] 0.755
[1] 0.76
[1] 0.765
[1] 0.77
[1] 0.775
[1] 0.78
[1] 0.785
[1] 0.79
[1] 0.795
[1] 0.8
[1] 0.805
[1] 0.81
[1] 0.815
[1] 0.82
[1] 0.825
[1] 0.83
[1] 0.835
[1] 0.84
[1] 0.845
[1] 0.85
[1] 0.855
[1] 0.86
[1] 0.865
[1] 0.87
[1] 0.875
[1] 0.88
[1] 0.885
[1] 0.89
[1] 0.895
[1] 0.9
[1] 0.905
[1] 0.91
[1] 0.915
[1] 0.92
[1] 0.925
[1] 0.93
[1] 0.935
[1] 0.94
[1] 0.945
[1] 0.95
[1] 0.955
[1] 0.96
[1] 0.965
[1] 0.97
[1] 0.975
[1] 0.98
[1] 0.985
[1] 0.99
[1] 0.995
[1] 1
if (sum(abs(Y) != abs(Re(Y))) == 0) {
  Y <- Re(Y)
}
nrow(input_var)
[1] 200
nrow(Y)
[1] 200
Y <- as.data.frame(Y); colnames(Y) <- c("eigenvalue2", "A_dens", "trace")
(data <- cbind(input_var, Y)) %>% head()

Plot Results

pairs(data)

cor(data)
                     p        beta
p            1.0000000  0.00000000
beta         0.0000000  1.00000000
size         0.0000000  0.00000000
eigenvalue2 -0.3657134  0.57594018
A_dens       0.9944206 -0.01047333
trace       -0.3986436 -0.70810602
                   size eigenvalue2
p            0.00000000  -0.3657134
beta         0.00000000   0.5759402
size         1.00000000  -0.3528876
eigenvalue2 -0.35288757   1.0000000
A_dens       0.05079861  -0.4078047
trace       -0.04276000  -0.4008003
                 A_dens      trace
p            0.99442058 -0.3986436
beta        -0.01047333 -0.7081060
size         0.05079861 -0.0427600
eigenvalue2 -0.40780467 -0.4008003
A_dens       1.00000000 -0.3935818
trace       -0.39358178  1.0000000
library(corrplot)
cormat = cor(data, method = 'spearman')
corrplot(cormat, method = "ellipse", type = "lower")

names(data)
[1] "p"           "beta"        "size"       
[4] "eigenvalue2" "A_dens"      "trace"      

Let’s look at a 3d Output

# plot3d(data$beta, data$A_dens, data$eigenvalue2, size = 2)

names(data)
[1] "p"           "beta"        "size"       
[4] "eigenvalue2" "A_dens"      "trace"      
library(plotly)

d <- data
# d$beta <- log(d$beta)

fig <- plot_ly(d, x = ~A_dens, y = ~beta, z = ~eigenvalue2)
fig <- fig %>% add_markers(size = 1)
fig <- fig %>% layout(scene = list(xaxis = list(title = 'Density'),
                     yaxis = list(title = 'Beta'),
                     zaxis = list(title = 'E2')))

fig

NA

names(data)
[1] "p"           "beta"        "size"       
[4] "eigenvalue2" "A_dens"      "trace"      
library(plotly)

#d <- data[sample(1:nrow(data), 1000),]
d <- data

fig <- plot_ly(d, x = ~A_dens, y = ~trace, z = ~eigenvalue2)
fig <- fig %>% add_markers(size = 1)
fig <- fig %>% layout(scene = list(xaxis = list(title = 'Density'),
                     yaxis = list(title = 'trace'),
                     zaxis = list(title = 'E2')))

fig

NA

Clearly I should be able to model this.

Let’s look at density, it’s continuous so I should be able to take splices of A_dens.

A_dens is dependent on p, so I’ll just use different ps

p    <- seq(from = 0.1, to = 0.95, length.out = 5)
beta <- seq(from = 1, to = 10, length.out = 1000)
size = 100
input_var <- expand.grid("p" = p, "beta" = beta, "size" = size)
input_var

nc <- length(random_graph(1, 1, 1))
Y <- matrix(ncol = nc, nrow = nrow(input_var))
for (i in 1:nrow(input_var)) {
  X <- as.vector(input_var[i,])
  Y[i,] <-  random_graph(X$p, X$beta, X$size)
  print(i/nrow(input_var))
}
[1] 2e-04
[1] 4e-04
[1] 6e-04
[1] 8e-04
[1] 0.001
[1] 0.0012
[1] 0.0014
[1] 0.0016
[1] 0.0018
[1] 0.002
[1] 0.0022
[1] 0.0024
[1] 0.0026
[1] 0.0028
[1] 0.003
[1] 0.0032
[1] 0.0034
[1] 0.0036
[1] 0.0038
[1] 0.004
[1] 0.0042
[1] 0.0044
[1] 0.0046
[1] 0.0048
[1] 0.005
[1] 0.0052
[1] 0.0054
[1] 0.0056
[1] 0.0058
[1] 0.006
[1] 0.0062
[1] 0.0064
[1] 0.0066
[1] 0.0068
[1] 0.007
[1] 0.0072
[1] 0.0074
[1] 0.0076
[1] 0.0078
[1] 0.008
[1] 0.0082
[1] 0.0084
[1] 0.0086
[1] 0.0088
[1] 0.009
[1] 0.0092
[1] 0.0094
[1] 0.0096
[1] 0.0098
[1] 0.01
[1] 0.0102
[1] 0.0104
[1] 0.0106
[1] 0.0108
[1] 0.011
[1] 0.0112
[1] 0.0114
[1] 0.0116
[1] 0.0118
[1] 0.012
[1] 0.0122
[1] 0.0124
[1] 0.0126
[1] 0.0128
[1] 0.013
[1] 0.0132
[1] 0.0134
[1] 0.0136
[1] 0.0138
[1] 0.014
[1] 0.0142
[1] 0.0144
[1] 0.0146
[1] 0.0148
[1] 0.015
[1] 0.0152
[1] 0.0154
[1] 0.0156
[1] 0.0158
[1] 0.016
[1] 0.0162
[1] 0.0164
[1] 0.0166
[1] 0.0168
[1] 0.017
[1] 0.0172
[1] 0.0174
[1] 0.0176
[1] 0.0178
[1] 0.018
[1] 0.0182
[1] 0.0184
[1] 0.0186
[1] 0.0188
[1] 0.019
[1] 0.0192
[1] 0.0194
[1] 0.0196
[1] 0.0198
[1] 0.02
[1] 0.0202
[1] 0.0204
[1] 0.0206
[1] 0.0208
[1] 0.021
[1] 0.0212
[1] 0.0214
[1] 0.0216
[1] 0.0218
[1] 0.022
[1] 0.0222
[1] 0.0224
[1] 0.0226
[1] 0.0228
[1] 0.023
[1] 0.0232
[1] 0.0234
[1] 0.0236
[1] 0.0238
[1] 0.024
[1] 0.0242
[1] 0.0244
[1] 0.0246
[1] 0.0248
[1] 0.025
[1] 0.0252
[1] 0.0254
[1] 0.0256
[1] 0.0258
[1] 0.026
[1] 0.0262
[1] 0.0264
[1] 0.0266
[1] 0.0268
[1] 0.027
[1] 0.0272
[1] 0.0274
[1] 0.0276
[1] 0.0278
[1] 0.028
[1] 0.0282
[1] 0.0284
[1] 0.0286
[1] 0.0288
[1] 0.029
[1] 0.0292
[1] 0.0294
[1] 0.0296
[1] 0.0298
[1] 0.03
[1] 0.0302
[1] 0.0304
[1] 0.0306
[1] 0.0308
[1] 0.031
[1] 0.0312
[1] 0.0314
[1] 0.0316
[1] 0.0318
[1] 0.032
[1] 0.0322
[1] 0.0324
[1] 0.0326
[1] 0.0328
[1] 0.033
[1] 0.0332
[1] 0.0334
[1] 0.0336
[1] 0.0338
[1] 0.034
[1] 0.0342
[1] 0.0344
[1] 0.0346
[1] 0.0348
[1] 0.035
[1] 0.0352
[1] 0.0354
[1] 0.0356
[1] 0.0358
[1] 0.036
[1] 0.0362
[1] 0.0364
[1] 0.0366
[1] 0.0368
[1] 0.037
[1] 0.0372
[1] 0.0374
[1] 0.0376
[1] 0.0378
[1] 0.038
[1] 0.0382
[1] 0.0384
[1] 0.0386
[1] 0.0388
[1] 0.039
[1] 0.0392
[1] 0.0394
[1] 0.0396
[1] 0.0398
[1] 0.04
[1] 0.0402
[1] 0.0404
[1] 0.0406
[1] 0.0408
[1] 0.041
[1] 0.0412
[1] 0.0414
[1] 0.0416
[1] 0.0418
[1] 0.042
[1] 0.0422
[1] 0.0424
[1] 0.0426
[1] 0.0428
[1] 0.043
[1] 0.0432
[1] 0.0434
[1] 0.0436
[1] 0.0438
[1] 0.044
[1] 0.0442
[1] 0.0444
[1] 0.0446
[1] 0.0448
[1] 0.045
[1] 0.0452
[1] 0.0454
[1] 0.0456
[1] 0.0458
[1] 0.046
[1] 0.0462
[1] 0.0464
[1] 0.0466
[1] 0.0468
[1] 0.047
[1] 0.0472
[1] 0.0474
[1] 0.0476
[1] 0.0478
[1] 0.048
[1] 0.0482
[1] 0.0484
[1] 0.0486
[1] 0.0488
[1] 0.049
[1] 0.0492
[1] 0.0494
[1] 0.0496
[1] 0.0498
[1] 0.05
[1] 0.0502
[1] 0.0504
[1] 0.0506
[1] 0.0508
[1] 0.051
[1] 0.0512
[1] 0.0514
[1] 0.0516
[1] 0.0518
[1] 0.052
[1] 0.0522
[1] 0.0524
[1] 0.0526
[1] 0.0528
[1] 0.053
[1] 0.0532
[1] 0.0534
[1] 0.0536
[1] 0.0538
[1] 0.054
[1] 0.0542
[1] 0.0544
[1] 0.0546
[1] 0.0548
[1] 0.055
[1] 0.0552
[1] 0.0554
[1] 0.0556
[1] 0.0558
[1] 0.056
[1] 0.0562
[1] 0.0564
[1] 0.0566
[1] 0.0568
[1] 0.057
[1] 0.0572
[1] 0.0574
[1] 0.0576
[1] 0.0578
[1] 0.058
[1] 0.0582
[1] 0.0584
[1] 0.0586
[1] 0.0588
[1] 0.059
[1] 0.0592
[1] 0.0594
[1] 0.0596
[1] 0.0598
[1] 0.06
[1] 0.0602
[1] 0.0604
[1] 0.0606
[1] 0.0608
[1] 0.061
[1] 0.0612
[1] 0.0614
[1] 0.0616
[1] 0.0618
[1] 0.062
[1] 0.0622
[1] 0.0624
[1] 0.0626
[1] 0.0628
[1] 0.063
[1] 0.0632
[1] 0.0634
[1] 0.0636
[1] 0.0638
[1] 0.064
[1] 0.0642
[1] 0.0644
[1] 0.0646
[1] 0.0648
[1] 0.065
[1] 0.0652
[1] 0.0654
[1] 0.0656
[1] 0.0658
[1] 0.066
[1] 0.0662
[1] 0.0664
[1] 0.0666
[1] 0.0668
[1] 0.067
[1] 0.0672
[1] 0.0674
[1] 0.0676
[1] 0.0678
[1] 0.068
[1] 0.0682
[1] 0.0684
[1] 0.0686
[1] 0.0688
[1] 0.069
[1] 0.0692
[1] 0.0694
[1] 0.0696
[1] 0.0698
[1] 0.07
[1] 0.0702
[1] 0.0704
[1] 0.0706
[1] 0.0708
[1] 0.071
[1] 0.0712
[1] 0.0714
[1] 0.0716
[1] 0.0718
[1] 0.072
[1] 0.0722
[1] 0.0724
[1] 0.0726
[1] 0.0728
[1] 0.073
[1] 0.0732
[1] 0.0734
[1] 0.0736
[1] 0.0738
[1] 0.074
[1] 0.0742
[1] 0.0744
[1] 0.0746
[1] 0.0748
[1] 0.075
[1] 0.0752
[1] 0.0754
[1] 0.0756
[1] 0.0758
[1] 0.076
[1] 0.0762
[1] 0.0764
[1] 0.0766
[1] 0.0768
[1] 0.077
[1] 0.0772
[1] 0.0774
[1] 0.0776
[1] 0.0778
[1] 0.078
[1] 0.0782
[1] 0.0784
[1] 0.0786
[1] 0.0788
[1] 0.079
[1] 0.0792
[1] 0.0794
[1] 0.0796
[1] 0.0798
[1] 0.08
[1] 0.0802
[1] 0.0804
[1] 0.0806
[1] 0.0808
[1] 0.081
[1] 0.0812
[1] 0.0814
[1] 0.0816
[1] 0.0818
[1] 0.082
[1] 0.0822
[1] 0.0824
[1] 0.0826
[1] 0.0828
[1] 0.083
[1] 0.0832
[1] 0.0834
[1] 0.0836
[1] 0.0838
[1] 0.084
[1] 0.0842
[1] 0.0844
[1] 0.0846
[1] 0.0848
[1] 0.085
[1] 0.0852
[1] 0.0854
[1] 0.0856
[1] 0.0858
[1] 0.086
[1] 0.0862
[1] 0.0864
[1] 0.0866
[1] 0.0868
[1] 0.087
[1] 0.0872
[1] 0.0874
[1] 0.0876
[1] 0.0878
[1] 0.088
[1] 0.0882
[1] 0.0884
[1] 0.0886
[1] 0.0888
[1] 0.089
[1] 0.0892
[1] 0.0894
[1] 0.0896
[1] 0.0898
[1] 0.09
[1] 0.0902
[1] 0.0904
[1] 0.0906
[1] 0.0908
[1] 0.091
[1] 0.0912
[1] 0.0914
[1] 0.0916
[1] 0.0918
[1] 0.092
[1] 0.0922
[1] 0.0924
[1] 0.0926
[1] 0.0928
[1] 0.093
[1] 0.0932
[1] 0.0934
[1] 0.0936
[1] 0.0938
[1] 0.094
[1] 0.0942
[1] 0.0944
[1] 0.0946
[1] 0.0948
[1] 0.095
[1] 0.0952
[1] 0.0954
[1] 0.0956
[1] 0.0958
[1] 0.096
[1] 0.0962
[1] 0.0964
[1] 0.0966
[1] 0.0968
[1] 0.097
[1] 0.0972
[1] 0.0974
[1] 0.0976
[1] 0.0978
[1] 0.098
[1] 0.0982
[1] 0.0984
[1] 0.0986
[1] 0.0988
[1] 0.099
[1] 0.0992
[1] 0.0994
[1] 0.0996
[1] 0.0998
[1] 0.1
[1] 0.1002
[1] 0.1004
[1] 0.1006
[1] 0.1008
[1] 0.101
[1] 0.1012
[1] 0.1014
[1] 0.1016
[1] 0.1018
[1] 0.102
[1] 0.1022
[1] 0.1024
[1] 0.1026
[1] 0.1028
[1] 0.103
[1] 0.1032
[1] 0.1034
[1] 0.1036
[1] 0.1038
[1] 0.104
[1] 0.1042
[1] 0.1044
[1] 0.1046
[1] 0.1048
[1] 0.105
[1] 0.1052
[1] 0.1054
[1] 0.1056
[1] 0.1058
[1] 0.106
[1] 0.1062
[1] 0.1064
[1] 0.1066
[1] 0.1068
[1] 0.107
[1] 0.1072
[1] 0.1074
[1] 0.1076
[1] 0.1078
[1] 0.108
[1] 0.1082
[1] 0.1084
[1] 0.1086
[1] 0.1088
[1] 0.109
[1] 0.1092
[1] 0.1094
[1] 0.1096
[1] 0.1098
[1] 0.11
[1] 0.1102
[1] 0.1104
[1] 0.1106
[1] 0.1108
[1] 0.111
[1] 0.1112
[1] 0.1114
[1] 0.1116
[1] 0.1118
[1] 0.112
[1] 0.1122
[1] 0.1124
[1] 0.1126
[1] 0.1128
[1] 0.113
[1] 0.1132
[1] 0.1134
[1] 0.1136
[1] 0.1138
[1] 0.114
[1] 0.1142
[1] 0.1144
[1] 0.1146
[1] 0.1148
[1] 0.115
[1] 0.1152
[1] 0.1154
[1] 0.1156
[1] 0.1158
[1] 0.116
[1] 0.1162
[1] 0.1164
[1] 0.1166
[1] 0.1168
[1] 0.117
[1] 0.1172
[1] 0.1174
[1] 0.1176
[1] 0.1178
[1] 0.118
[1] 0.1182
[1] 0.1184
[1] 0.1186
[1] 0.1188
[1] 0.119
[1] 0.1192
[1] 0.1194
[1] 0.1196
[1] 0.1198
[1] 0.12
[1] 0.1202
[1] 0.1204
[1] 0.1206
[1] 0.1208
[1] 0.121
[1] 0.1212
[1] 0.1214
[1] 0.1216
[1] 0.1218
[1] 0.122
[1] 0.1222
[1] 0.1224
[1] 0.1226
[1] 0.1228
[1] 0.123
[1] 0.1232
[1] 0.1234
[1] 0.1236
[1] 0.1238
[1] 0.124
[1] 0.1242
[1] 0.1244
[1] 0.1246
[1] 0.1248
[1] 0.125
[1] 0.1252
[1] 0.1254
[1] 0.1256
[1] 0.1258
[1] 0.126
[1] 0.1262
[1] 0.1264
[1] 0.1266
[1] 0.1268
[1] 0.127
[1] 0.1272
[1] 0.1274
[1] 0.1276
[1] 0.1278
[1] 0.128
[1] 0.1282
[1] 0.1284
[1] 0.1286
[1] 0.1288
[1] 0.129
[1] 0.1292
[1] 0.1294
[1] 0.1296
[1] 0.1298
[1] 0.13
[1] 0.1302
[1] 0.1304
[1] 0.1306
[1] 0.1308
[1] 0.131
[1] 0.1312
[1] 0.1314
[1] 0.1316
[1] 0.1318
[1] 0.132
[1] 0.1322
[1] 0.1324
[1] 0.1326
[1] 0.1328
[1] 0.133
[1] 0.1332
[1] 0.1334
[1] 0.1336
[1] 0.1338
[1] 0.134
[1] 0.1342
[1] 0.1344
[1] 0.1346
[1] 0.1348
[1] 0.135
[1] 0.1352
[1] 0.1354
[1] 0.1356
[1] 0.1358
[1] 0.136
[1] 0.1362
[1] 0.1364
[1] 0.1366
[1] 0.1368
[1] 0.137
[1] 0.1372
[1] 0.1374
[1] 0.1376
[1] 0.1378
[1] 0.138
[1] 0.1382
[1] 0.1384
[1] 0.1386
[1] 0.1388
[1] 0.139
[1] 0.1392
[1] 0.1394
[1] 0.1396
[1] 0.1398
[1] 0.14
[1] 0.1402
[1] 0.1404
[1] 0.1406
[1] 0.1408
[1] 0.141
[1] 0.1412
[1] 0.1414
[1] 0.1416
[1] 0.1418
[1] 0.142
[1] 0.1422
[1] 0.1424
[1] 0.1426
[1] 0.1428
[1] 0.143
[1] 0.1432
[1] 0.1434
[1] 0.1436
[1] 0.1438
[1] 0.144
[1] 0.1442
[1] 0.1444
[1] 0.1446
[1] 0.1448
[1] 0.145
[1] 0.1452
[1] 0.1454
[1] 0.1456
[1] 0.1458
[1] 0.146
[1] 0.1462
[1] 0.1464
[1] 0.1466
[1] 0.1468
[1] 0.147
[1] 0.1472
[1] 0.1474
[1] 0.1476
[1] 0.1478
[1] 0.148
[1] 0.1482
[1] 0.1484
[1] 0.1486
[1] 0.1488
[1] 0.149
[1] 0.1492
[1] 0.1494
[1] 0.1496
[1] 0.1498
[1] 0.15
[1] 0.1502
[1] 0.1504
[1] 0.1506
[1] 0.1508
[1] 0.151
[1] 0.1512
[1] 0.1514
[1] 0.1516
[1] 0.1518
[1] 0.152
[1] 0.1522
[1] 0.1524
[1] 0.1526
[1] 0.1528
[1] 0.153
[1] 0.1532
[1] 0.1534
[1] 0.1536
[1] 0.1538
[1] 0.154
[1] 0.1542
[1] 0.1544
[1] 0.1546
[1] 0.1548
[1] 0.155
[1] 0.1552
[1] 0.1554
[1] 0.1556
[1] 0.1558
[1] 0.156
[1] 0.1562
[1] 0.1564
[1] 0.1566
[1] 0.1568
[1] 0.157
[1] 0.1572
[1] 0.1574
[1] 0.1576
[1] 0.1578
[1] 0.158
[1] 0.1582
[1] 0.1584
[1] 0.1586
[1] 0.1588
[1] 0.159
[1] 0.1592
[1] 0.1594
[1] 0.1596
[1] 0.1598
[1] 0.16
[1] 0.1602
[1] 0.1604
[1] 0.1606
[1] 0.1608
[1] 0.161
[1] 0.1612
[1] 0.1614
[1] 0.1616
[1] 0.1618
[1] 0.162
[1] 0.1622
[1] 0.1624
[1] 0.1626
[1] 0.1628
[1] 0.163
[1] 0.1632
[1] 0.1634
[1] 0.1636
[1] 0.1638
[1] 0.164
[1] 0.1642
[1] 0.1644
[1] 0.1646
[1] 0.1648
[1] 0.165
[1] 0.1652
[1] 0.1654
[1] 0.1656
[1] 0.1658
[1] 0.166
[1] 0.1662
[1] 0.1664
[1] 0.1666
[1] 0.1668
[1] 0.167
[1] 0.1672
[1] 0.1674
[1] 0.1676
[1] 0.1678
[1] 0.168
[1] 0.1682
[1] 0.1684
[1] 0.1686
[1] 0.1688
[1] 0.169
[1] 0.1692
[1] 0.1694
[1] 0.1696
[1] 0.1698
[1] 0.17
[1] 0.1702
[1] 0.1704
[1] 0.1706
[1] 0.1708
[1] 0.171
[1] 0.1712
[1] 0.1714
[1] 0.1716
[1] 0.1718
[1] 0.172
[1] 0.1722
[1] 0.1724
[1] 0.1726
[1] 0.1728
[1] 0.173
[1] 0.1732
[1] 0.1734
[1] 0.1736
[1] 0.1738
[1] 0.174
[1] 0.1742
[1] 0.1744
[1] 0.1746
[1] 0.1748
[1] 0.175
[1] 0.1752
[1] 0.1754
[1] 0.1756
[1] 0.1758
[1] 0.176
[1] 0.1762
[1] 0.1764
[1] 0.1766
[1] 0.1768
[1] 0.177
[1] 0.1772
[1] 0.1774
[1] 0.1776
[1] 0.1778
[1] 0.178
[1] 0.1782
[1] 0.1784
[1] 0.1786
[1] 0.1788
[1] 0.179
[1] 0.1792
[1] 0.1794
[1] 0.1796
[1] 0.1798
[1] 0.18
[1] 0.1802
[1] 0.1804
[1] 0.1806
[1] 0.1808
[1] 0.181
[1] 0.1812
[1] 0.1814
[1] 0.1816
[1] 0.1818
[1] 0.182
[1] 0.1822
[1] 0.1824
[1] 0.1826
[1] 0.1828
[1] 0.183
[1] 0.1832
[1] 0.1834
[1] 0.1836
[1] 0.1838
[1] 0.184
[1] 0.1842
[1] 0.1844
[1] 0.1846
[1] 0.1848
[1] 0.185
[1] 0.1852
[1] 0.1854
[1] 0.1856
[1] 0.1858
[1] 0.186
[1] 0.1862
[1] 0.1864
[1] 0.1866
[1] 0.1868
[1] 0.187
[1] 0.1872
[1] 0.1874
[1] 0.1876
[1] 0.1878
[1] 0.188
[1] 0.1882
[1] 0.1884
[1] 0.1886
[1] 0.1888
[1] 0.189
[1] 0.1892
[1] 0.1894
[1] 0.1896
[1] 0.1898
[1] 0.19
[1] 0.1902
[1] 0.1904
[1] 0.1906
[1] 0.1908
[1] 0.191
[1] 0.1912
[1] 0.1914
[1] 0.1916
[1] 0.1918
[1] 0.192
[1] 0.1922
[1] 0.1924
[1] 0.1926
[1] 0.1928
[1] 0.193
[1] 0.1932
[1] 0.1934
[1] 0.1936
[1] 0.1938
[1] 0.194
[1] 0.1942
[1] 0.1944
[1] 0.1946
[1] 0.1948
[1] 0.195
[1] 0.1952
[1] 0.1954
[1] 0.1956
[1] 0.1958
[1] 0.196
[1] 0.1962
[1] 0.1964
[1] 0.1966
[1] 0.1968
[1] 0.197
[1] 0.1972
[1] 0.1974
[1] 0.1976
[1] 0.1978
[1] 0.198
[1] 0.1982
[1] 0.1984
[1] 0.1986
[1] 0.1988
[1] 0.199
[1] 0.1992
[1] 0.1994
[1] 0.1996
[1] 0.1998
[1] 0.2
[1] 0.2002
[1] 0.2004
[1] 0.2006
[1] 0.2008
[1] 0.201
[1] 0.2012
[1] 0.2014
[1] 0.2016
[1] 0.2018
[1] 0.202
[1] 0.2022
[1] 0.2024
[1] 0.2026
[1] 0.2028
[1] 0.203
[1] 0.2032
[1] 0.2034
[1] 0.2036
[1] 0.2038
[1] 0.204
[1] 0.2042
[1] 0.2044
[1] 0.2046
[1] 0.2048
[1] 0.205
[1] 0.2052
[1] 0.2054
[1] 0.2056
[1] 0.2058
[1] 0.206
[1] 0.2062
[1] 0.2064
[1] 0.2066
[1] 0.2068
[1] 0.207
[1] 0.2072
[1] 0.2074
[1] 0.2076
[1] 0.2078
[1] 0.208
[1] 0.2082
[1] 0.2084
[1] 0.2086
[1] 0.2088
[1] 0.209
[1] 0.2092
[1] 0.2094
[1] 0.2096
[1] 0.2098
[1] 0.21
[1] 0.2102
[1] 0.2104
[1] 0.2106
[1] 0.2108
[1] 0.211
[1] 0.2112
[1] 0.2114
[1] 0.2116
[1] 0.2118
[1] 0.212
[1] 0.2122
[1] 0.2124
[1] 0.2126
[1] 0.2128
[1] 0.213
[1] 0.2132
[1] 0.2134
[1] 0.2136
[1] 0.2138
[1] 0.214
[1] 0.2142
[1] 0.2144
[1] 0.2146
[1] 0.2148
[1] 0.215
[1] 0.2152
[1] 0.2154
[1] 0.2156
[1] 0.2158
[1] 0.216
[1] 0.2162
[1] 0.2164
[1] 0.2166
[1] 0.2168
[1] 0.217
[1] 0.2172
[1] 0.2174
[1] 0.2176
[1] 0.2178
[1] 0.218
[1] 0.2182
[1] 0.2184
[1] 0.2186
[1] 0.2188
[1] 0.219
[1] 0.2192
[1] 0.2194
[1] 0.2196
[1] 0.2198
[1] 0.22
[1] 0.2202
[1] 0.2204
[1] 0.2206
[1] 0.2208
[1] 0.221
[1] 0.2212
[1] 0.2214
[1] 0.2216
[1] 0.2218
[1] 0.222
[1] 0.2222
[1] 0.2224
[1] 0.2226
[1] 0.2228
[1] 0.223
[1] 0.2232
[1] 0.2234
[1] 0.2236
[1] 0.2238
[1] 0.224
[1] 0.2242
[1] 0.2244
[1] 0.2246
[1] 0.2248
[1] 0.225
[1] 0.2252
[1] 0.2254
[1] 0.2256
[1] 0.2258
[1] 0.226
[1] 0.2262
[1] 0.2264
[1] 0.2266
[1] 0.2268
[1] 0.227
[1] 0.2272
[1] 0.2274
[1] 0.2276
[1] 0.2278
[1] 0.228
[1] 0.2282
[1] 0.2284
[1] 0.2286
[1] 0.2288
[1] 0.229
[1] 0.2292
[1] 0.2294
[1] 0.2296
[1] 0.2298
[1] 0.23
[1] 0.2302
[1] 0.2304
[1] 0.2306
[1] 0.2308
[1] 0.231
[1] 0.2312
[1] 0.2314
[1] 0.2316
[1] 0.2318
[1] 0.232
[1] 0.2322
[1] 0.2324
[1] 0.2326
[1] 0.2328
[1] 0.233
[1] 0.2332
[1] 0.2334
[1] 0.2336
[1] 0.2338
[1] 0.234
[1] 0.2342
[1] 0.2344
[1] 0.2346
[1] 0.2348
[1] 0.235
[1] 0.2352
[1] 0.2354
[1] 0.2356
[1] 0.2358
[1] 0.236
[1] 0.2362
[1] 0.2364
[1] 0.2366
[1] 0.2368
[1] 0.237
[1] 0.2372
[1] 0.2374
[1] 0.2376
[1] 0.2378
[1] 0.238
[1] 0.2382
[1] 0.2384
[1] 0.2386
[1] 0.2388
[1] 0.239
[1] 0.2392
[1] 0.2394
[1] 0.2396
[1] 0.2398
[1] 0.24
[1] 0.2402
[1] 0.2404
[1] 0.2406
[1] 0.2408
[1] 0.241
[1] 0.2412
[1] 0.2414
[1] 0.2416
[1] 0.2418
[1] 0.242
[1] 0.2422
[1] 0.2424
[1] 0.2426
[1] 0.2428
[1] 0.243
[1] 0.2432
[1] 0.2434
[1] 0.2436
[1] 0.2438
[1] 0.244
[1] 0.2442
[1] 0.2444
[1] 0.2446
[1] 0.2448
[1] 0.245
[1] 0.2452
[1] 0.2454
[1] 0.2456
[1] 0.2458
[1] 0.246
[1] 0.2462
[1] 0.2464
[1] 0.2466
[1] 0.2468
[1] 0.247
[1] 0.2472
[1] 0.2474
[1] 0.2476
[1] 0.2478
[1] 0.248
[1] 0.2482
[1] 0.2484
[1] 0.2486
[1] 0.2488
[1] 0.249
[1] 0.2492
[1] 0.2494
[1] 0.2496
[1] 0.2498
[1] 0.25
[1] 0.2502
[1] 0.2504
[1] 0.2506
[1] 0.2508
[1] 0.251
[1] 0.2512
[1] 0.2514
[1] 0.2516
[1] 0.2518
[1] 0.252
[1] 0.2522
[1] 0.2524
[1] 0.2526
[1] 0.2528
[1] 0.253
[1] 0.2532
[1] 0.2534
[1] 0.2536
[1] 0.2538
[1] 0.254
[1] 0.2542
[1] 0.2544
[1] 0.2546
[1] 0.2548
[1] 0.255
[1] 0.2552
[1] 0.2554
[1] 0.2556
[1] 0.2558
[1] 0.256
[1] 0.2562
[1] 0.2564
[1] 0.2566
[1] 0.2568
[1] 0.257
[1] 0.2572
[1] 0.2574
[1] 0.2576
[1] 0.2578
[1] 0.258
[1] 0.2582
[1] 0.2584
[1] 0.2586
[1] 0.2588
[1] 0.259
[1] 0.2592
[1] 0.2594
[1] 0.2596
[1] 0.2598
[1] 0.26
[1] 0.2602
[1] 0.2604
[1] 0.2606
[1] 0.2608
[1] 0.261
[1] 0.2612
[1] 0.2614
[1] 0.2616
[1] 0.2618
[1] 0.262
[1] 0.2622
[1] 0.2624
[1] 0.2626
[1] 0.2628
[1] 0.263
[1] 0.2632
[1] 0.2634
[1] 0.2636
[1] 0.2638
[1] 0.264
[1] 0.2642
[1] 0.2644
[1] 0.2646
[1] 0.2648
[1] 0.265
[1] 0.2652
[1] 0.2654
[1] 0.2656
[1] 0.2658
[1] 0.266
[1] 0.2662
[1] 0.2664
[1] 0.2666
[1] 0.2668
[1] 0.267
[1] 0.2672
[1] 0.2674
[1] 0.2676
[1] 0.2678
[1] 0.268
[1] 0.2682
[1] 0.2684
[1] 0.2686
[1] 0.2688
[1] 0.269
[1] 0.2692
[1] 0.2694
[1] 0.2696
[1] 0.2698
[1] 0.27
[1] 0.2702
[1] 0.2704
[1] 0.2706
[1] 0.2708
[1] 0.271
[1] 0.2712
[1] 0.2714
[1] 0.2716
[1] 0.2718
[1] 0.272
[1] 0.2722
[1] 0.2724
[1] 0.2726
[1] 0.2728
[1] 0.273
[1] 0.2732
[1] 0.2734
[1] 0.2736
[1] 0.2738
[1] 0.274
[1] 0.2742
[1] 0.2744
[1] 0.2746
[1] 0.2748
[1] 0.275
[1] 0.2752
[1] 0.2754
[1] 0.2756
[1] 0.2758
[1] 0.276
[1] 0.2762
[1] 0.2764
[1] 0.2766
[1] 0.2768
[1] 0.277
[1] 0.2772
[1] 0.2774
[1] 0.2776
[1] 0.2778
[1] 0.278
[1] 0.2782
[1] 0.2784
[1] 0.2786
[1] 0.2788
[1] 0.279
[1] 0.2792
[1] 0.2794
[1] 0.2796
[1] 0.2798
[1] 0.28
[1] 0.2802
[1] 0.2804
[1] 0.2806
[1] 0.2808
[1] 0.281
[1] 0.2812
[1] 0.2814
[1] 0.2816
[1] 0.2818
[1] 0.282
[1] 0.2822
[1] 0.2824
[1] 0.2826
[1] 0.2828
[1] 0.283
[1] 0.2832
[1] 0.2834
[1] 0.2836
[1] 0.2838
[1] 0.284
[1] 0.2842
[1] 0.2844
[1] 0.2846
[1] 0.2848
[1] 0.285
[1] 0.2852
[1] 0.2854
[1] 0.2856
[1] 0.2858
[1] 0.286
[1] 0.2862
[1] 0.2864
[1] 0.2866
[1] 0.2868
[1] 0.287
[1] 0.2872
[1] 0.2874
[1] 0.2876
[1] 0.2878
[1] 0.288
[1] 0.2882
[1] 0.2884
[1] 0.2886
[1] 0.2888
[1] 0.289
[1] 0.2892
[1] 0.2894
[1] 0.2896
[1] 0.2898
[1] 0.29
[1] 0.2902
[1] 0.2904
[1] 0.2906
[1] 0.2908
[1] 0.291
[1] 0.2912
[1] 0.2914
[1] 0.2916
[1] 0.2918
[1] 0.292
[1] 0.2922
[1] 0.2924
[1] 0.2926
[1] 0.2928
[1] 0.293
[1] 0.2932
[1] 0.2934
[1] 0.2936
[1] 0.2938
[1] 0.294
[1] 0.2942
[1] 0.2944
[1] 0.2946
[1] 0.2948
[1] 0.295
[1] 0.2952
[1] 0.2954
[1] 0.2956
[1] 0.2958
[1] 0.296
[1] 0.2962
[1] 0.2964
[1] 0.2966
[1] 0.2968
[1] 0.297
[1] 0.2972
[1] 0.2974
[1] 0.2976
[1] 0.2978
[1] 0.298
[1] 0.2982
[1] 0.2984
[1] 0.2986
[1] 0.2988
[1] 0.299
[1] 0.2992
[1] 0.2994
[1] 0.2996
[1] 0.2998
[1] 0.3
[1] 0.3002
[1] 0.3004
[1] 0.3006
[1] 0.3008
[1] 0.301
[1] 0.3012
[1] 0.3014
[1] 0.3016
[1] 0.3018
[1] 0.302
[1] 0.3022
[1] 0.3024
[1] 0.3026
[1] 0.3028
[1] 0.303
[1] 0.3032
[1] 0.3034
[1] 0.3036
[1] 0.3038
[1] 0.304
[1] 0.3042
[1] 0.3044
[1] 0.3046
[1] 0.3048
[1] 0.305
[1] 0.3052
[1] 0.3054
[1] 0.3056
[1] 0.3058
[1] 0.306
[1] 0.3062
[1] 0.3064
[1] 0.3066
[1] 0.3068
[1] 0.307
[1] 0.3072
[1] 0.3074
[1] 0.3076
[1] 0.3078
[1] 0.308
[1] 0.3082
[1] 0.3084
[1] 0.3086
[1] 0.3088
[1] 0.309
[1] 0.3092
[1] 0.3094
[1] 0.3096
[1] 0.3098
[1] 0.31
[1] 0.3102
[1] 0.3104
[1] 0.3106
[1] 0.3108
[1] 0.311
[1] 0.3112
[1] 0.3114
[1] 0.3116
[1] 0.3118
[1] 0.312
[1] 0.3122
[1] 0.3124
[1] 0.3126
[1] 0.3128
[1] 0.313
[1] 0.3132
[1] 0.3134
[1] 0.3136
[1] 0.3138
[1] 0.314
[1] 0.3142
[1] 0.3144
[1] 0.3146
[1] 0.3148
[1] 0.315
[1] 0.3152
[1] 0.3154
[1] 0.3156
[1] 0.3158
[1] 0.316
[1] 0.3162
[1] 0.3164
[1] 0.3166
[1] 0.3168
[1] 0.317
[1] 0.3172
[1] 0.3174
[1] 0.3176
[1] 0.3178
[1] 0.318
[1] 0.3182
[1] 0.3184
[1] 0.3186
[1] 0.3188
[1] 0.319
[1] 0.3192
[1] 0.3194
[1] 0.3196
[1] 0.3198
[1] 0.32
[1] 0.3202
[1] 0.3204
[1] 0.3206
[1] 0.3208
[1] 0.321
[1] 0.3212
[1] 0.3214
[1] 0.3216
[1] 0.3218
[1] 0.322
[1] 0.3222
[1] 0.3224
[1] 0.3226
[1] 0.3228
[1] 0.323
[1] 0.3232
[1] 0.3234
[1] 0.3236
[1] 0.3238
[1] 0.324
[1] 0.3242
[1] 0.3244
[1] 0.3246
[1] 0.3248
[1] 0.325
[1] 0.3252
[1] 0.3254
[1] 0.3256
[1] 0.3258
[1] 0.326
[1] 0.3262
[1] 0.3264
[1] 0.3266
[1] 0.3268
[1] 0.327
[1] 0.3272
[1] 0.3274
[1] 0.3276
[1] 0.3278
[1] 0.328
[1] 0.3282
[1] 0.3284
[1] 0.3286
[1] 0.3288
[1] 0.329
[1] 0.3292
[1] 0.3294
[1] 0.3296
[1] 0.3298
[1] 0.33
[1] 0.3302
[1] 0.3304
[1] 0.3306
[1] 0.3308
[1] 0.331
[1] 0.3312
[1] 0.3314
[1] 0.3316
[1] 0.3318
[1] 0.332
[1] 0.3322
[1] 0.3324
[1] 0.3326
[1] 0.3328
[1] 0.333
[1] 0.3332
[1] 0.3334
[1] 0.3336
[1] 0.3338
[1] 0.334
[1] 0.3342
[1] 0.3344
[1] 0.3346
[1] 0.3348
[1] 0.335
[1] 0.3352
[1] 0.3354
[1] 0.3356
[1] 0.3358
[1] 0.336
[1] 0.3362
[1] 0.3364
[1] 0.3366
[1] 0.3368
[1] 0.337
[1] 0.3372
[1] 0.3374
[1] 0.3376
[1] 0.3378
[1] 0.338
[1] 0.3382
[1] 0.3384
[1] 0.3386
[1] 0.3388
[1] 0.339
[1] 0.3392
[1] 0.3394
[1] 0.3396
[1] 0.3398
[1] 0.34
[1] 0.3402
[1] 0.3404
[1] 0.3406
[1] 0.3408
[1] 0.341
[1] 0.3412
[1] 0.3414
[1] 0.3416
[1] 0.3418
[1] 0.342
[1] 0.3422
[1] 0.3424
[1] 0.3426
[1] 0.3428
[1] 0.343
[1] 0.3432
[1] 0.3434
[1] 0.3436
[1] 0.3438
[1] 0.344
[1] 0.3442
[1] 0.3444
[1] 0.3446
[1] 0.3448
[1] 0.345
[1] 0.3452
[1] 0.3454
[1] 0.3456
[1] 0.3458
[1] 0.346
[1] 0.3462
[1] 0.3464
[1] 0.3466
[1] 0.3468
[1] 0.347
[1] 0.3472
[1] 0.3474
[1] 0.3476
[1] 0.3478
[1] 0.348
[1] 0.3482
[1] 0.3484
[1] 0.3486
[1] 0.3488
[1] 0.349
[1] 0.3492
[1] 0.3494
[1] 0.3496
[1] 0.3498
[1] 0.35
[1] 0.3502
[1] 0.3504
[1] 0.3506
[1] 0.3508
[1] 0.351
[1] 0.3512
[1] 0.3514
[1] 0.3516
[1] 0.3518
[1] 0.352
[1] 0.3522
[1] 0.3524
[1] 0.3526
[1] 0.3528
[1] 0.353
[1] 0.3532
[1] 0.3534
[1] 0.3536
[1] 0.3538
[1] 0.354
[1] 0.3542
[1] 0.3544
[1] 0.3546
[1] 0.3548
[1] 0.355
[1] 0.3552
[1] 0.3554
[1] 0.3556
[1] 0.3558
[1] 0.356
[1] 0.3562
[1] 0.3564
[1] 0.3566
[1] 0.3568
[1] 0.357
[1] 0.3572
[1] 0.3574
[1] 0.3576
[1] 0.3578
[1] 0.358
[1] 0.3582
[1] 0.3584
[1] 0.3586
[1] 0.3588
[1] 0.359
[1] 0.3592
[1] 0.3594
[1] 0.3596
[1] 0.3598
[1] 0.36
[1] 0.3602
[1] 0.3604
[1] 0.3606
[1] 0.3608
[1] 0.361
[1] 0.3612
[1] 0.3614
[1] 0.3616
[1] 0.3618
[1] 0.362
[1] 0.3622
[1] 0.3624
[1] 0.3626
[1] 0.3628
[1] 0.363
[1] 0.3632
[1] 0.3634
[1] 0.3636
[1] 0.3638
[1] 0.364
[1] 0.3642
[1] 0.3644
[1] 0.3646
[1] 0.3648
[1] 0.365
[1] 0.3652
[1] 0.3654
[1] 0.3656
[1] 0.3658
[1] 0.366
[1] 0.3662
[1] 0.3664
[1] 0.3666
[1] 0.3668
[1] 0.367
[1] 0.3672
[1] 0.3674
[1] 0.3676
[1] 0.3678
[1] 0.368
[1] 0.3682
[1] 0.3684
[1] 0.3686
[1] 0.3688
[1] 0.369
[1] 0.3692
[1] 0.3694
[1] 0.3696
[1] 0.3698
[1] 0.37
[1] 0.3702
[1] 0.3704
[1] 0.3706
[1] 0.3708
[1] 0.371
[1] 0.3712
[1] 0.3714
[1] 0.3716
[1] 0.3718
[1] 0.372
[1] 0.3722
[1] 0.3724
[1] 0.3726
[1] 0.3728
[1] 0.373
[1] 0.3732
[1] 0.3734
[1] 0.3736
[1] 0.3738
[1] 0.374
[1] 0.3742
[1] 0.3744
[1] 0.3746
[1] 0.3748
[1] 0.375
[1] 0.3752
[1] 0.3754
[1] 0.3756
[1] 0.3758
[1] 0.376
[1] 0.3762
[1] 0.3764
[1] 0.3766
[1] 0.3768
[1] 0.377
[1] 0.3772
[1] 0.3774
[1] 0.3776
[1] 0.3778
[1] 0.378
[1] 0.3782
[1] 0.3784
[1] 0.3786
[1] 0.3788
[1] 0.379
[1] 0.3792
[1] 0.3794
[1] 0.3796
[1] 0.3798
[1] 0.38
[1] 0.3802
[1] 0.3804
[1] 0.3806
[1] 0.3808
[1] 0.381
[1] 0.3812
[1] 0.3814
[1] 0.3816
[1] 0.3818
[1] 0.382
[1] 0.3822
[1] 0.3824
[1] 0.3826
[1] 0.3828
[1] 0.383
[1] 0.3832
[1] 0.3834
[1] 0.3836
[1] 0.3838
[1] 0.384
[1] 0.3842
[1] 0.3844
[1] 0.3846
[1] 0.3848
[1] 0.385
[1] 0.3852
[1] 0.3854
[1] 0.3856
[1] 0.3858
[1] 0.386
[1] 0.3862
[1] 0.3864
[1] 0.3866
[1] 0.3868
[1] 0.387
[1] 0.3872
[1] 0.3874
[1] 0.3876
[1] 0.3878
[1] 0.388
[1] 0.3882
[1] 0.3884
[1] 0.3886
[1] 0.3888
[1] 0.389
[1] 0.3892
[1] 0.3894
[1] 0.3896
[1] 0.3898
[1] 0.39
[1] 0.3902
[1] 0.3904
[1] 0.3906
[1] 0.3908
[1] 0.391
[1] 0.3912
[1] 0.3914
[1] 0.3916
[1] 0.3918
[1] 0.392
[1] 0.3922
[1] 0.3924
[1] 0.3926
[1] 0.3928
[1] 0.393
[1] 0.3932
[1] 0.3934
[1] 0.3936
[1] 0.3938
[1] 0.394
[1] 0.3942
[1] 0.3944
[1] 0.3946
[1] 0.3948
[1] 0.395
[1] 0.3952
[1] 0.3954
[1] 0.3956
[1] 0.3958
[1] 0.396
[1] 0.3962
[1] 0.3964
[1] 0.3966
[1] 0.3968
[1] 0.397
[1] 0.3972
[1] 0.3974
[1] 0.3976
[1] 0.3978
[1] 0.398
[1] 0.3982
[1] 0.3984
[1] 0.3986
[1] 0.3988
[1] 0.399
[1] 0.3992
[1] 0.3994
[1] 0.3996
[1] 0.3998
[1] 0.4
[1] 0.4002
[1] 0.4004
[1] 0.4006
[1] 0.4008
[1] 0.401
[1] 0.4012
[1] 0.4014
[1] 0.4016
[1] 0.4018
[1] 0.402
[1] 0.4022
[1] 0.4024
[1] 0.4026
[1] 0.4028
[1] 0.403
[1] 0.4032
[1] 0.4034
[1] 0.4036
[1] 0.4038
[1] 0.404
[1] 0.4042
[1] 0.4044
[1] 0.4046
[1] 0.4048
[1] 0.405
[1] 0.4052
[1] 0.4054
[1] 0.4056
[1] 0.4058
[1] 0.406
[1] 0.4062
[1] 0.4064
[1] 0.4066
[1] 0.4068
[1] 0.407
[1] 0.4072
[1] 0.4074
[1] 0.4076
[1] 0.4078
[1] 0.408
[1] 0.4082
[1] 0.4084
[1] 0.4086
[1] 0.4088
[1] 0.409
[1] 0.4092
[1] 0.4094
[1] 0.4096
[1] 0.4098
[1] 0.41
[1] 0.4102
[1] 0.4104
[1] 0.4106
[1] 0.4108
[1] 0.411
[1] 0.4112
[1] 0.4114
[1] 0.4116
[1] 0.4118
[1] 0.412
[1] 0.4122
[1] 0.4124
[1] 0.4126
[1] 0.4128
[1] 0.413
[1] 0.4132
[1] 0.4134
[1] 0.4136
[1] 0.4138
[1] 0.414
[1] 0.4142
[1] 0.4144
[1] 0.4146
[1] 0.4148
[1] 0.415
[1] 0.4152
[1] 0.4154
[1] 0.4156
[1] 0.4158
[1] 0.416
[1] 0.4162
[1] 0.4164
[1] 0.4166
[1] 0.4168
[1] 0.417
[1] 0.4172
[1] 0.4174
[1] 0.4176
[1] 0.4178
[1] 0.418
[1] 0.4182
[1] 0.4184
[1] 0.4186
[1] 0.4188
[1] 0.419
[1] 0.4192
[1] 0.4194
[1] 0.4196
[1] 0.4198
[1] 0.42
[1] 0.4202
[1] 0.4204
[1] 0.4206
[1] 0.4208
[1] 0.421
[1] 0.4212
[1] 0.4214
[1] 0.4216
[1] 0.4218
[1] 0.422
[1] 0.4222
[1] 0.4224
[1] 0.4226
[1] 0.4228
[1] 0.423
[1] 0.4232
[1] 0.4234
[1] 0.4236
[1] 0.4238
[1] 0.424
[1] 0.4242
[1] 0.4244
[1] 0.4246
[1] 0.4248
[1] 0.425
[1] 0.4252
[1] 0.4254
[1] 0.4256
[1] 0.4258
[1] 0.426
[1] 0.4262
[1] 0.4264
[1] 0.4266
[1] 0.4268
[1] 0.427
[1] 0.4272
[1] 0.4274
[1] 0.4276
[1] 0.4278
[1] 0.428
[1] 0.4282
[1] 0.4284
[1] 0.4286
[1] 0.4288
[1] 0.429
[1] 0.4292
[1] 0.4294
[1] 0.4296
[1] 0.4298
[1] 0.43
[1] 0.4302
[1] 0.4304
[1] 0.4306
[1] 0.4308
[1] 0.431
[1] 0.4312
[1] 0.4314
[1] 0.4316
[1] 0.4318
[1] 0.432
[1] 0.4322
[1] 0.4324
[1] 0.4326
[1] 0.4328
[1] 0.433
[1] 0.4332
[1] 0.4334
[1] 0.4336
[1] 0.4338
[1] 0.434
[1] 0.4342
[1] 0.4344
[1] 0.4346
[1] 0.4348
[1] 0.435
[1] 0.4352
[1] 0.4354
[1] 0.4356
[1] 0.4358
[1] 0.436
[1] 0.4362
[1] 0.4364
[1] 0.4366
[1] 0.4368
[1] 0.437
[1] 0.4372
[1] 0.4374
[1] 0.4376
[1] 0.4378
[1] 0.438
[1] 0.4382
[1] 0.4384
[1] 0.4386
[1] 0.4388
[1] 0.439
[1] 0.4392
[1] 0.4394
[1] 0.4396
[1] 0.4398
[1] 0.44
[1] 0.4402
[1] 0.4404
[1] 0.4406
[1] 0.4408
[1] 0.441
[1] 0.4412
[1] 0.4414
[1] 0.4416
[1] 0.4418
[1] 0.442
[1] 0.4422
[1] 0.4424
[1] 0.4426
[1] 0.4428
[1] 0.443
[1] 0.4432
[1] 0.4434
[1] 0.4436
[1] 0.4438
[1] 0.444
[1] 0.4442
[1] 0.4444
[1] 0.4446
[1] 0.4448
[1] 0.445
[1] 0.4452
[1] 0.4454
[1] 0.4456
[1] 0.4458
[1] 0.446
[1] 0.4462
[1] 0.4464
[1] 0.4466
[1] 0.4468
[1] 0.447
[1] 0.4472
[1] 0.4474
[1] 0.4476
[1] 0.4478
[1] 0.448
[1] 0.4482
[1] 0.4484
[1] 0.4486
[1] 0.4488
[1] 0.449
[1] 0.4492
[1] 0.4494
[1] 0.4496
[1] 0.4498
[1] 0.45
[1] 0.4502
[1] 0.4504
[1] 0.4506
[1] 0.4508
[1] 0.451
[1] 0.4512
[1] 0.4514
[1] 0.4516
[1] 0.4518
[1] 0.452
[1] 0.4522
[1] 0.4524
[1] 0.4526
[1] 0.4528
[1] 0.453
[1] 0.4532
[1] 0.4534
[1] 0.4536
[1] 0.4538
[1] 0.454
[1] 0.4542
[1] 0.4544
[1] 0.4546
[1] 0.4548
[1] 0.455
[1] 0.4552
[1] 0.4554
[1] 0.4556
[1] 0.4558
[1] 0.456
[1] 0.4562
[1] 0.4564
[1] 0.4566
[1] 0.4568
[1] 0.457
[1] 0.4572
[1] 0.4574
[1] 0.4576
[1] 0.4578
[1] 0.458
[1] 0.4582
[1] 0.4584
[1] 0.4586
[1] 0.4588
[1] 0.459
[1] 0.4592
[1] 0.4594
[1] 0.4596
[1] 0.4598
[1] 0.46
[1] 0.4602
[1] 0.4604
[1] 0.4606
[1] 0.4608
[1] 0.461
[1] 0.4612
[1] 0.4614
[1] 0.4616
[1] 0.4618
[1] 0.462
[1] 0.4622
[1] 0.4624
[1] 0.4626
[1] 0.4628
[1] 0.463
[1] 0.4632
[1] 0.4634
[1] 0.4636
[1] 0.4638
[1] 0.464
[1] 0.4642
[1] 0.4644
[1] 0.4646
[1] 0.4648
[1] 0.465
[1] 0.4652
[1] 0.4654
[1] 0.4656
[1] 0.4658
[1] 0.466
[1] 0.4662
[1] 0.4664
[1] 0.4666
[1] 0.4668
[1] 0.467
[1] 0.4672
[1] 0.4674
[1] 0.4676
[1] 0.4678
[1] 0.468
[1] 0.4682
[1] 0.4684
[1] 0.4686
[1] 0.4688
[1] 0.469
[1] 0.4692
[1] 0.4694
[1] 0.4696
[1] 0.4698
[1] 0.47
[1] 0.4702
[1] 0.4704
[1] 0.4706
[1] 0.4708
[1] 0.471
[1] 0.4712
[1] 0.4714
[1] 0.4716
[1] 0.4718
[1] 0.472
[1] 0.4722
[1] 0.4724
[1] 0.4726
[1] 0.4728
[1] 0.473
[1] 0.4732
[1] 0.4734
[1] 0.4736
[1] 0.4738
[1] 0.474
[1] 0.4742
[1] 0.4744
[1] 0.4746
[1] 0.4748
[1] 0.475
[1] 0.4752
[1] 0.4754
[1] 0.4756
[1] 0.4758
[1] 0.476
[1] 0.4762
[1] 0.4764
[1] 0.4766
[1] 0.4768
[1] 0.477
[1] 0.4772
[1] 0.4774
[1] 0.4776
[1] 0.4778
[1] 0.478
[1] 0.4782
[1] 0.4784
[1] 0.4786
[1] 0.4788
[1] 0.479
[1] 0.4792
[1] 0.4794
[1] 0.4796
[1] 0.4798
[1] 0.48
[1] 0.4802
[1] 0.4804
[1] 0.4806
[1] 0.4808
[1] 0.481
[1] 0.4812
[1] 0.4814
[1] 0.4816
[1] 0.4818
[1] 0.482
[1] 0.4822
[1] 0.4824
[1] 0.4826
[1] 0.4828
[1] 0.483
[1] 0.4832
[1] 0.4834
[1] 0.4836
[1] 0.4838
[1] 0.484
[1] 0.4842
[1] 0.4844
[1] 0.4846
[1] 0.4848
[1] 0.485
[1] 0.4852
[1] 0.4854
[1] 0.4856
[1] 0.4858
[1] 0.486
[1] 0.4862
[1] 0.4864
[1] 0.4866
[1] 0.4868
[1] 0.487
[1] 0.4872
[1] 0.4874
[1] 0.4876
[1] 0.4878
[1] 0.488
[1] 0.4882
[1] 0.4884
[1] 0.4886
[1] 0.4888
[1] 0.489
[1] 0.4892
[1] 0.4894
[1] 0.4896
[1] 0.4898
[1] 0.49
[1] 0.4902
[1] 0.4904
[1] 0.4906
[1] 0.4908
[1] 0.491
[1] 0.4912
[1] 0.4914
[1] 0.4916
[1] 0.4918
[1] 0.492
[1] 0.4922
[1] 0.4924
[1] 0.4926
[1] 0.4928
[1] 0.493
[1] 0.4932
[1] 0.4934
[1] 0.4936
[1] 0.4938
[1] 0.494
[1] 0.4942
[1] 0.4944
[1] 0.4946
[1] 0.4948
[1] 0.495
[1] 0.4952
[1] 0.4954
[1] 0.4956
[1] 0.4958
[1] 0.496
[1] 0.4962
[1] 0.4964
[1] 0.4966
[1] 0.4968
[1] 0.497
[1] 0.4972
[1] 0.4974
[1] 0.4976
[1] 0.4978
[1] 0.498
[1] 0.4982
[1] 0.4984
[1] 0.4986
[1] 0.4988
[1] 0.499
[1] 0.4992
[1] 0.4994
[1] 0.4996
[1] 0.4998
[1] 0.5
[1] 0.5002
[1] 0.5004
[1] 0.5006
[1] 0.5008
[1] 0.501
[1] 0.5012
[1] 0.5014
[1] 0.5016
[1] 0.5018
[1] 0.502
[1] 0.5022
[1] 0.5024
[1] 0.5026
[1] 0.5028
[1] 0.503
[1] 0.5032
[1] 0.5034
[1] 0.5036
[1] 0.5038
[1] 0.504
[1] 0.5042
[1] 0.5044
[1] 0.5046
[1] 0.5048
[1] 0.505
[1] 0.5052
[1] 0.5054
[1] 0.5056
[1] 0.5058
[1] 0.506
[1] 0.5062
[1] 0.5064
[1] 0.5066
[1] 0.5068
[1] 0.507
[1] 0.5072
[1] 0.5074
[1] 0.5076
[1] 0.5078
[1] 0.508
[1] 0.5082
[1] 0.5084
[1] 0.5086
[1] 0.5088
[1] 0.509
[1] 0.5092
[1] 0.5094
[1] 0.5096
[1] 0.5098
[1] 0.51
[1] 0.5102
[1] 0.5104
[1] 0.5106
[1] 0.5108
[1] 0.511
[1] 0.5112
[1] 0.5114
[1] 0.5116
[1] 0.5118
[1] 0.512
[1] 0.5122
[1] 0.5124
[1] 0.5126
[1] 0.5128
[1] 0.513
[1] 0.5132
[1] 0.5134
[1] 0.5136
[1] 0.5138
[1] 0.514
[1] 0.5142
[1] 0.5144
[1] 0.5146
[1] 0.5148
[1] 0.515
[1] 0.5152
[1] 0.5154
[1] 0.5156
[1] 0.5158
[1] 0.516
[1] 0.5162
[1] 0.5164
[1] 0.5166
[1] 0.5168
[1] 0.517
[1] 0.5172
[1] 0.5174
[1] 0.5176
[1] 0.5178
[1] 0.518
[1] 0.5182
[1] 0.5184
[1] 0.5186
[1] 0.5188
[1] 0.519
[1] 0.5192
[1] 0.5194
[1] 0.5196
[1] 0.5198
[1] 0.52
[1] 0.5202
[1] 0.5204
[1] 0.5206
[1] 0.5208
[1] 0.521
[1] 0.5212
[1] 0.5214
[1] 0.5216
[1] 0.5218
[1] 0.522
[1] 0.5222
[1] 0.5224
[1] 0.5226
[1] 0.5228
[1] 0.523
[1] 0.5232
[1] 0.5234
[1] 0.5236
[1] 0.5238
[1] 0.524
[1] 0.5242
[1] 0.5244
[1] 0.5246
[1] 0.5248
[1] 0.525
[1] 0.5252
[1] 0.5254
[1] 0.5256
[1] 0.5258
[1] 0.526
[1] 0.5262
[1] 0.5264
[1] 0.5266
[1] 0.5268
[1] 0.527
[1] 0.5272
[1] 0.5274
[1] 0.5276
[1] 0.5278
[1] 0.528
[1] 0.5282
[1] 0.5284
[1] 0.5286
[1] 0.5288
[1] 0.529
[1] 0.5292
[1] 0.5294
[1] 0.5296
[1] 0.5298
[1] 0.53
[1] 0.5302
[1] 0.5304
[1] 0.5306
[1] 0.5308
[1] 0.531
[1] 0.5312
[1] 0.5314
[1] 0.5316
[1] 0.5318
[1] 0.532
[1] 0.5322
[1] 0.5324
[1] 0.5326
[1] 0.5328
[1] 0.533
[1] 0.5332
[1] 0.5334
[1] 0.5336
[1] 0.5338
[1] 0.534
[1] 0.5342
[1] 0.5344
[1] 0.5346
[1] 0.5348
[1] 0.535
[1] 0.5352
[1] 0.5354
[1] 0.5356
[1] 0.5358
[1] 0.536
[1] 0.5362
[1] 0.5364
[1] 0.5366
[1] 0.5368
[1] 0.537
[1] 0.5372
[1] 0.5374
[1] 0.5376
[1] 0.5378
[1] 0.538
[1] 0.5382
[1] 0.5384
[1] 0.5386
[1] 0.5388
[1] 0.539
[1] 0.5392
[1] 0.5394
[1] 0.5396
[1] 0.5398
[1] 0.54
[1] 0.5402
[1] 0.5404
[1] 0.5406
[1] 0.5408
[1] 0.541
[1] 0.5412
[1] 0.5414
[1] 0.5416
[1] 0.5418
[1] 0.542
[1] 0.5422
[1] 0.5424
[1] 0.5426
[1] 0.5428
[1] 0.543
[1] 0.5432
[1] 0.5434
[1] 0.5436
[1] 0.5438
[1] 0.544
[1] 0.5442
[1] 0.5444
[1] 0.5446
[1] 0.5448
[1] 0.545
[1] 0.5452
[1] 0.5454
[1] 0.5456
[1] 0.5458
[1] 0.546
[1] 0.5462
[1] 0.5464
[1] 0.5466
[1] 0.5468
[1] 0.547
[1] 0.5472
[1] 0.5474
[1] 0.5476
[1] 0.5478
[1] 0.548
[1] 0.5482
[1] 0.5484
[1] 0.5486
[1] 0.5488
[1] 0.549
[1] 0.5492
[1] 0.5494
[1] 0.5496
[1] 0.5498
[1] 0.55
[1] 0.5502
[1] 0.5504
[1] 0.5506
[1] 0.5508
[1] 0.551
[1] 0.5512
[1] 0.5514
[1] 0.5516
[1] 0.5518
[1] 0.552
[1] 0.5522
[1] 0.5524
[1] 0.5526
[1] 0.5528
[1] 0.553
[1] 0.5532
[1] 0.5534
[1] 0.5536
[1] 0.5538
[1] 0.554
[1] 0.5542
[1] 0.5544
[1] 0.5546
[1] 0.5548
[1] 0.555
[1] 0.5552
[1] 0.5554
[1] 0.5556
[1] 0.5558
[1] 0.556
[1] 0.5562
[1] 0.5564
[1] 0.5566
[1] 0.5568
[1] 0.557
[1] 0.5572
[1] 0.5574
[1] 0.5576
[1] 0.5578
[1] 0.558
[1] 0.5582
[1] 0.5584
[1] 0.5586
[1] 0.5588
[1] 0.559
[1] 0.5592
[1] 0.5594
[1] 0.5596
[1] 0.5598
[1] 0.56
[1] 0.5602
[1] 0.5604
[1] 0.5606
[1] 0.5608
[1] 0.561
[1] 0.5612
[1] 0.5614
[1] 0.5616
[1] 0.5618
[1] 0.562
[1] 0.5622
[1] 0.5624
[1] 0.5626
[1] 0.5628
[1] 0.563
[1] 0.5632
[1] 0.5634
[1] 0.5636
[1] 0.5638
[1] 0.564
[1] 0.5642
[1] 0.5644
[1] 0.5646
[1] 0.5648
[1] 0.565
[1] 0.5652
[1] 0.5654
[1] 0.5656
[1] 0.5658
[1] 0.566
[1] 0.5662
[1] 0.5664
[1] 0.5666
[1] 0.5668
[1] 0.567
[1] 0.5672
[1] 0.5674
[1] 0.5676
[1] 0.5678
[1] 0.568
[1] 0.5682
[1] 0.5684
[1] 0.5686
[1] 0.5688
[1] 0.569
[1] 0.5692
[1] 0.5694
[1] 0.5696
[1] 0.5698
[1] 0.57
[1] 0.5702
[1] 0.5704
[1] 0.5706
[1] 0.5708
[1] 0.571
[1] 0.5712
[1] 0.5714
[1] 0.5716
[1] 0.5718
[1] 0.572
[1] 0.5722
[1] 0.5724
[1] 0.5726
[1] 0.5728
[1] 0.573
[1] 0.5732
[1] 0.5734
[1] 0.5736
[1] 0.5738
[1] 0.574
[1] 0.5742
[1] 0.5744
[1] 0.5746
[1] 0.5748
[1] 0.575
[1] 0.5752
[1] 0.5754
[1] 0.5756
[1] 0.5758
[1] 0.576
[1] 0.5762
[1] 0.5764
[1] 0.5766
[1] 0.5768
[1] 0.577
[1] 0.5772
[1] 0.5774
[1] 0.5776
[1] 0.5778
[1] 0.578
[1] 0.5782
[1] 0.5784
[1] 0.5786
[1] 0.5788
[1] 0.579
[1] 0.5792
[1] 0.5794
[1] 0.5796
[1] 0.5798
[1] 0.58
[1] 0.5802
[1] 0.5804
[1] 0.5806
[1] 0.5808
[1] 0.581
[1] 0.5812
[1] 0.5814
[1] 0.5816
[1] 0.5818
[1] 0.582
[1] 0.5822
[1] 0.5824
[1] 0.5826
[1] 0.5828
[1] 0.583
[1] 0.5832
[1] 0.5834
[1] 0.5836
[1] 0.5838
[1] 0.584
[1] 0.5842
[1] 0.5844
[1] 0.5846
[1] 0.5848
[1] 0.585
[1] 0.5852
[1] 0.5854
[1] 0.5856
[1] 0.5858
[1] 0.586
[1] 0.5862
[1] 0.5864
[1] 0.5866
[1] 0.5868
[1] 0.587
[1] 0.5872
[1] 0.5874
[1] 0.5876
[1] 0.5878
[1] 0.588
[1] 0.5882
[1] 0.5884
[1] 0.5886
[1] 0.5888
[1] 0.589
[1] 0.5892
[1] 0.5894
[1] 0.5896
[1] 0.5898
[1] 0.59
[1] 0.5902
[1] 0.5904
[1] 0.5906
[1] 0.5908
[1] 0.591
[1] 0.5912
[1] 0.5914
[1] 0.5916
[1] 0.5918
[1] 0.592
[1] 0.5922
[1] 0.5924
[1] 0.5926
[1] 0.5928
[1] 0.593
[1] 0.5932
[1] 0.5934
[1] 0.5936
[1] 0.5938
[1] 0.594
[1] 0.5942
[1] 0.5944
[1] 0.5946
[1] 0.5948
[1] 0.595
[1] 0.5952
[1] 0.5954
[1] 0.5956
[1] 0.5958
[1] 0.596
[1] 0.5962
[1] 0.5964
[1] 0.5966
[1] 0.5968
[1] 0.597
[1] 0.5972
[1] 0.5974
[1] 0.5976
[1] 0.5978
[1] 0.598
[1] 0.5982
[1] 0.5984
[1] 0.5986
[1] 0.5988
[1] 0.599
[1] 0.5992
[1] 0.5994
[1] 0.5996
[1] 0.5998
[1] 0.6
[1] 0.6002
[1] 0.6004
[1] 0.6006
[1] 0.6008
[1] 0.601
[1] 0.6012
[1] 0.6014
[1] 0.6016
[1] 0.6018
[1] 0.602
[1] 0.6022
[1] 0.6024
[1] 0.6026
[1] 0.6028
[1] 0.603
[1] 0.6032
[1] 0.6034
[1] 0.6036
[1] 0.6038
[1] 0.604
[1] 0.6042
[1] 0.6044
[1] 0.6046
[1] 0.6048
[1] 0.605
[1] 0.6052
[1] 0.6054
[1] 0.6056
[1] 0.6058
[1] 0.606
[1] 0.6062
[1] 0.6064
[1] 0.6066
[1] 0.6068
[1] 0.607
[1] 0.6072
[1] 0.6074
[1] 0.6076
[1] 0.6078
[1] 0.608
[1] 0.6082
[1] 0.6084
[1] 0.6086
[1] 0.6088
[1] 0.609
[1] 0.6092
[1] 0.6094
[1] 0.6096
[1] 0.6098
[1] 0.61
[1] 0.6102
[1] 0.6104
[1] 0.6106
[1] 0.6108
[1] 0.611
[1] 0.6112
[1] 0.6114
[1] 0.6116
[1] 0.6118
[1] 0.612
[1] 0.6122
[1] 0.6124
[1] 0.6126
[1] 0.6128
[1] 0.613
[1] 0.6132
[1] 0.6134
[1] 0.6136
[1] 0.6138
[1] 0.614
[1] 0.6142
[1] 0.6144
[1] 0.6146
[1] 0.6148
[1] 0.615
[1] 0.6152
[1] 0.6154
[1] 0.6156
[1] 0.6158
[1] 0.616
[1] 0.6162
[1] 0.6164
[1] 0.6166
[1] 0.6168
[1] 0.617
[1] 0.6172
[1] 0.6174
[1] 0.6176
[1] 0.6178
[1] 0.618
[1] 0.6182
[1] 0.6184
[1] 0.6186
[1] 0.6188
[1] 0.619
[1] 0.6192
[1] 0.6194
[1] 0.6196
[1] 0.6198
[1] 0.62
[1] 0.6202
[1] 0.6204
[1] 0.6206
[1] 0.6208
[1] 0.621
[1] 0.6212
[1] 0.6214
[1] 0.6216
[1] 0.6218
[1] 0.622
[1] 0.6222
[1] 0.6224
[1] 0.6226
[1] 0.6228
[1] 0.623
[1] 0.6232
[1] 0.6234
[1] 0.6236
[1] 0.6238
[1] 0.624
[1] 0.6242
[1] 0.6244
[1] 0.6246
[1] 0.6248
[1] 0.625
[1] 0.6252
[1] 0.6254
[1] 0.6256
[1] 0.6258
[1] 0.626
[1] 0.6262
[1] 0.6264
[1] 0.6266
[1] 0.6268
[1] 0.627
[1] 0.6272
[1] 0.6274
[1] 0.6276
[1] 0.6278
[1] 0.628
[1] 0.6282
[1] 0.6284
[1] 0.6286
[1] 0.6288
[1] 0.629
[1] 0.6292
[1] 0.6294
[1] 0.6296
[1] 0.6298
[1] 0.63
[1] 0.6302
[1] 0.6304
[1] 0.6306
[1] 0.6308
[1] 0.631
[1] 0.6312
[1] 0.6314
[1] 0.6316
[1] 0.6318
[1] 0.632
[1] 0.6322
[1] 0.6324
[1] 0.6326
[1] 0.6328
[1] 0.633
[1] 0.6332
[1] 0.6334
[1] 0.6336
[1] 0.6338
[1] 0.634
[1] 0.6342
[1] 0.6344
[1] 0.6346
[1] 0.6348
[1] 0.635
[1] 0.6352
[1] 0.6354
[1] 0.6356
[1] 0.6358
[1] 0.636
[1] 0.6362
[1] 0.6364
[1] 0.6366
[1] 0.6368
[1] 0.637
[1] 0.6372
[1] 0.6374
[1] 0.6376
[1] 0.6378
[1] 0.638
[1] 0.6382
[1] 0.6384
[1] 0.6386
[1] 0.6388
[1] 0.639
[1] 0.6392
[1] 0.6394
[1] 0.6396
[1] 0.6398
[1] 0.64
[1] 0.6402
[1] 0.6404
[1] 0.6406
[1] 0.6408
[1] 0.641
[1] 0.6412
[1] 0.6414
[1] 0.6416
[1] 0.6418
[1] 0.642
[1] 0.6422
[1] 0.6424
[1] 0.6426
[1] 0.6428
[1] 0.643
[1] 0.6432
[1] 0.6434
[1] 0.6436
[1] 0.6438
[1] 0.644
[1] 0.6442
[1] 0.6444
[1] 0.6446
[1] 0.6448
[1] 0.645
[1] 0.6452
[1] 0.6454
[1] 0.6456
[1] 0.6458
[1] 0.646
[1] 0.6462
[1] 0.6464
[1] 0.6466
[1] 0.6468
[1] 0.647
[1] 0.6472
[1] 0.6474
[1] 0.6476
[1] 0.6478
[1] 0.648
[1] 0.6482
[1] 0.6484
[1] 0.6486
[1] 0.6488
[1] 0.649
[1] 0.6492
[1] 0.6494
[1] 0.6496
[1] 0.6498
[1] 0.65
[1] 0.6502
[1] 0.6504
[1] 0.6506
[1] 0.6508
[1] 0.651
[1] 0.6512
[1] 0.6514
[1] 0.6516
[1] 0.6518
[1] 0.652
[1] 0.6522
[1] 0.6524
[1] 0.6526
[1] 0.6528
[1] 0.653
[1] 0.6532
[1] 0.6534
[1] 0.6536
[1] 0.6538
[1] 0.654
[1] 0.6542
[1] 0.6544
[1] 0.6546
[1] 0.6548
[1] 0.655
[1] 0.6552
[1] 0.6554
[1] 0.6556
[1] 0.6558
[1] 0.656
[1] 0.6562
[1] 0.6564
[1] 0.6566
[1] 0.6568
[1] 0.657
[1] 0.6572
[1] 0.6574
[1] 0.6576
[1] 0.6578
[1] 0.658
[1] 0.6582
[1] 0.6584
[1] 0.6586
[1] 0.6588
[1] 0.659
[1] 0.6592
[1] 0.6594
[1] 0.6596
[1] 0.6598
[1] 0.66
[1] 0.6602
[1] 0.6604
[1] 0.6606
[1] 0.6608
[1] 0.661
[1] 0.6612
[1] 0.6614
[1] 0.6616
[1] 0.6618
[1] 0.662
[1] 0.6622
[1] 0.6624
[1] 0.6626
[1] 0.6628
[1] 0.663
[1] 0.6632
[1] 0.6634
[1] 0.6636
[1] 0.6638
[1] 0.664
[1] 0.6642
[1] 0.6644
[1] 0.6646
[1] 0.6648
[1] 0.665
[1] 0.6652
[1] 0.6654
[1] 0.6656
[1] 0.6658
[1] 0.666
[1] 0.6662
[1] 0.6664
[1] 0.6666
[1] 0.6668
[1] 0.667
[1] 0.6672
[1] 0.6674
[1] 0.6676
[1] 0.6678
[1] 0.668
[1] 0.6682
[1] 0.6684
[1] 0.6686
[1] 0.6688
[1] 0.669
[1] 0.6692
[1] 0.6694
[1] 0.6696
[1] 0.6698
[1] 0.67
[1] 0.6702
[1] 0.6704
[1] 0.6706
[1] 0.6708
[1] 0.671
[1] 0.6712
[1] 0.6714
[1] 0.6716
[1] 0.6718
[1] 0.672
[1] 0.6722
[1] 0.6724
[1] 0.6726
[1] 0.6728
[1] 0.673
[1] 0.6732
[1] 0.6734
[1] 0.6736
[1] 0.6738
[1] 0.674
[1] 0.6742
[1] 0.6744
[1] 0.6746
[1] 0.6748
[1] 0.675
[1] 0.6752
[1] 0.6754
[1] 0.6756
[1] 0.6758
[1] 0.676
[1] 0.6762
[1] 0.6764
[1] 0.6766
[1] 0.6768
[1] 0.677
[1] 0.6772
[1] 0.6774
[1] 0.6776
[1] 0.6778
[1] 0.678
[1] 0.6782
[1] 0.6784
[1] 0.6786
[1] 0.6788
[1] 0.679
[1] 0.6792
[1] 0.6794
[1] 0.6796
[1] 0.6798
[1] 0.68
[1] 0.6802
[1] 0.6804
[1] 0.6806
[1] 0.6808
[1] 0.681
[1] 0.6812
[1] 0.6814
[1] 0.6816
[1] 0.6818
[1] 0.682
[1] 0.6822
[1] 0.6824
[1] 0.6826
[1] 0.6828
[1] 0.683
[1] 0.6832
[1] 0.6834
[1] 0.6836
[1] 0.6838
[1] 0.684
[1] 0.6842
[1] 0.6844
[1] 0.6846
[1] 0.6848
[1] 0.685
[1] 0.6852
[1] 0.6854
[1] 0.6856
[1] 0.6858
[1] 0.686
[1] 0.6862
[1] 0.6864
[1] 0.6866
[1] 0.6868
[1] 0.687
[1] 0.6872
[1] 0.6874
[1] 0.6876
[1] 0.6878
[1] 0.688
[1] 0.6882
[1] 0.6884
[1] 0.6886
[1] 0.6888
[1] 0.689
[1] 0.6892
[1] 0.6894
[1] 0.6896
[1] 0.6898
[1] 0.69
[1] 0.6902
[1] 0.6904
[1] 0.6906
[1] 0.6908
[1] 0.691
[1] 0.6912
[1] 0.6914
[1] 0.6916
[1] 0.6918
[1] 0.692
[1] 0.6922
[1] 0.6924
[1] 0.6926
[1] 0.6928
[1] 0.693
[1] 0.6932
[1] 0.6934
[1] 0.6936
[1] 0.6938
[1] 0.694
[1] 0.6942
[1] 0.6944
[1] 0.6946
[1] 0.6948
[1] 0.695
[1] 0.6952
[1] 0.6954
[1] 0.6956
[1] 0.6958
[1] 0.696
[1] 0.6962
[1] 0.6964
[1] 0.6966
[1] 0.6968
[1] 0.697
[1] 0.6972
[1] 0.6974
[1] 0.6976
[1] 0.6978
[1] 0.698
[1] 0.6982
[1] 0.6984
[1] 0.6986
[1] 0.6988
[1] 0.699
[1] 0.6992
[1] 0.6994
[1] 0.6996
[1] 0.6998
[1] 0.7
[1] 0.7002
[1] 0.7004
[1] 0.7006
[1] 0.7008
[1] 0.701
[1] 0.7012
[1] 0.7014
[1] 0.7016
[1] 0.7018
[1] 0.702
[1] 0.7022
[1] 0.7024
[1] 0.7026
[1] 0.7028
[1] 0.703
[1] 0.7032
[1] 0.7034
[1] 0.7036
[1] 0.7038
[1] 0.704
[1] 0.7042
[1] 0.7044
[1] 0.7046
[1] 0.7048
[1] 0.705
[1] 0.7052
[1] 0.7054
[1] 0.7056
[1] 0.7058
[1] 0.706
[1] 0.7062
[1] 0.7064
[1] 0.7066
[1] 0.7068
[1] 0.707
[1] 0.7072
[1] 0.7074
[1] 0.7076
[1] 0.7078
[1] 0.708
[1] 0.7082
[1] 0.7084
[1] 0.7086
[1] 0.7088
[1] 0.709
[1] 0.7092
[1] 0.7094
[1] 0.7096
[1] 0.7098
[1] 0.71
[1] 0.7102
[1] 0.7104
[1] 0.7106
[1] 0.7108
[1] 0.711
[1] 0.7112
[1] 0.7114
[1] 0.7116
[1] 0.7118
[1] 0.712
[1] 0.7122
[1] 0.7124
[1] 0.7126
[1] 0.7128
[1] 0.713
[1] 0.7132
[1] 0.7134
[1] 0.7136
[1] 0.7138
[1] 0.714
[1] 0.7142
[1] 0.7144
[1] 0.7146
[1] 0.7148
[1] 0.715
[1] 0.7152
[1] 0.7154
[1] 0.7156
[1] 0.7158
[1] 0.716
[1] 0.7162
[1] 0.7164
[1] 0.7166
[1] 0.7168
[1] 0.717
[1] 0.7172
[1] 0.7174
[1] 0.7176
[1] 0.7178
[1] 0.718
[1] 0.7182
[1] 0.7184
[1] 0.7186
[1] 0.7188
[1] 0.719
[1] 0.7192
[1] 0.7194
[1] 0.7196
[1] 0.7198
[1] 0.72
[1] 0.7202
[1] 0.7204
[1] 0.7206
[1] 0.7208
[1] 0.721
[1] 0.7212
[1] 0.7214
[1] 0.7216
[1] 0.7218
[1] 0.722
[1] 0.7222
[1] 0.7224
[1] 0.7226
[1] 0.7228
[1] 0.723
[1] 0.7232
[1] 0.7234
[1] 0.7236
[1] 0.7238
[1] 0.724
[1] 0.7242
[1] 0.7244
[1] 0.7246
[1] 0.7248
[1] 0.725
[1] 0.7252
[1] 0.7254
[1] 0.7256
[1] 0.7258
[1] 0.726
[1] 0.7262
[1] 0.7264
[1] 0.7266
[1] 0.7268
[1] 0.727
[1] 0.7272
[1] 0.7274
[1] 0.7276
[1] 0.7278
[1] 0.728
[1] 0.7282
[1] 0.7284
[1] 0.7286
[1] 0.7288
[1] 0.729
[1] 0.7292
[1] 0.7294
[1] 0.7296
[1] 0.7298
[1] 0.73
[1] 0.7302
[1] 0.7304
[1] 0.7306
[1] 0.7308
[1] 0.731
[1] 0.7312
[1] 0.7314
[1] 0.7316
[1] 0.7318
[1] 0.732
[1] 0.7322
[1] 0.7324
[1] 0.7326
[1] 0.7328
[1] 0.733
[1] 0.7332
[1] 0.7334
[1] 0.7336
[1] 0.7338
[1] 0.734
[1] 0.7342
[1] 0.7344
[1] 0.7346
[1] 0.7348
[1] 0.735
[1] 0.7352
[1] 0.7354
[1] 0.7356
[1] 0.7358
[1] 0.736
[1] 0.7362
[1] 0.7364
[1] 0.7366
[1] 0.7368
[1] 0.737
[1] 0.7372
[1] 0.7374
[1] 0.7376
[1] 0.7378
[1] 0.738
[1] 0.7382
[1] 0.7384
[1] 0.7386
[1] 0.7388
[1] 0.739
[1] 0.7392
[1] 0.7394
[1] 0.7396
[1] 0.7398
[1] 0.74
[1] 0.7402
[1] 0.7404
[1] 0.7406
[1] 0.7408
[1] 0.741
[1] 0.7412
[1] 0.7414
[1] 0.7416
[1] 0.7418
[1] 0.742
[1] 0.7422
[1] 0.7424
[1] 0.7426
[1] 0.7428
[1] 0.743
[1] 0.7432
[1] 0.7434
[1] 0.7436
[1] 0.7438
[1] 0.744
[1] 0.7442
[1] 0.7444
[1] 0.7446
[1] 0.7448
[1] 0.745
[1] 0.7452
[1] 0.7454
[1] 0.7456
[1] 0.7458
[1] 0.746
[1] 0.7462
[1] 0.7464
[1] 0.7466
[1] 0.7468
[1] 0.747
[1] 0.7472
[1] 0.7474
[1] 0.7476
[1] 0.7478
[1] 0.748
[1] 0.7482
[1] 0.7484
[1] 0.7486
[1] 0.7488
[1] 0.749
[1] 0.7492
[1] 0.7494
[1] 0.7496
[1] 0.7498
[1] 0.75
[1] 0.7502
[1] 0.7504
[1] 0.7506
[1] 0.7508
[1] 0.751
[1] 0.7512
[1] 0.7514
[1] 0.7516
[1] 0.7518
[1] 0.752
[1] 0.7522
[1] 0.7524
[1] 0.7526
[1] 0.7528
[1] 0.753
[1] 0.7532
[1] 0.7534
[1] 0.7536
[1] 0.7538
[1] 0.754
[1] 0.7542
[1] 0.7544
[1] 0.7546
[1] 0.7548
[1] 0.755
[1] 0.7552
[1] 0.7554
[1] 0.7556
[1] 0.7558
[1] 0.756
[1] 0.7562
[1] 0.7564
[1] 0.7566
[1] 0.7568
[1] 0.757
[1] 0.7572
[1] 0.7574
[1] 0.7576
[1] 0.7578
[1] 0.758
[1] 0.7582
[1] 0.7584
[1] 0.7586
[1] 0.7588
[1] 0.759
[1] 0.7592
[1] 0.7594
[1] 0.7596
[1] 0.7598
[1] 0.76
[1] 0.7602
[1] 0.7604
[1] 0.7606
[1] 0.7608
[1] 0.761
[1] 0.7612
[1] 0.7614
[1] 0.7616
[1] 0.7618
[1] 0.762
[1] 0.7622
[1] 0.7624
[1] 0.7626
[1] 0.7628
[1] 0.763
[1] 0.7632
[1] 0.7634
[1] 0.7636
[1] 0.7638
[1] 0.764
[1] 0.7642
[1] 0.7644
[1] 0.7646
[1] 0.7648
[1] 0.765
[1] 0.7652
[1] 0.7654
[1] 0.7656
[1] 0.7658
[1] 0.766
[1] 0.7662
[1] 0.7664
[1] 0.7666
[1] 0.7668
[1] 0.767
[1] 0.7672
[1] 0.7674
[1] 0.7676
[1] 0.7678
[1] 0.768
[1] 0.7682
[1] 0.7684
[1] 0.7686
[1] 0.7688
[1] 0.769
[1] 0.7692
[1] 0.7694
[1] 0.7696
[1] 0.7698
[1] 0.77
[1] 0.7702
[1] 0.7704
[1] 0.7706
[1] 0.7708
[1] 0.771
[1] 0.7712
[1] 0.7714
[1] 0.7716
[1] 0.7718
[1] 0.772
[1] 0.7722
[1] 0.7724
[1] 0.7726
[1] 0.7728
[1] 0.773
[1] 0.7732
[1] 0.7734
[1] 0.7736
[1] 0.7738
[1] 0.774
[1] 0.7742
[1] 0.7744
[1] 0.7746
[1] 0.7748
[1] 0.775
[1] 0.7752
[1] 0.7754
[1] 0.7756
[1] 0.7758
[1] 0.776
[1] 0.7762
[1] 0.7764
[1] 0.7766
[1] 0.7768
[1] 0.777
[1] 0.7772
[1] 0.7774
[1] 0.7776
[1] 0.7778
[1] 0.778
[1] 0.7782
[1] 0.7784
[1] 0.7786
[1] 0.7788
[1] 0.779
[1] 0.7792
[1] 0.7794
[1] 0.7796
[1] 0.7798
[1] 0.78
[1] 0.7802
[1] 0.7804
[1] 0.7806
[1] 0.7808
[1] 0.781
[1] 0.7812
[1] 0.7814
[1] 0.7816
[1] 0.7818
[1] 0.782
[1] 0.7822
[1] 0.7824
[1] 0.7826
[1] 0.7828
[1] 0.783
[1] 0.7832
[1] 0.7834
[1] 0.7836
[1] 0.7838
[1] 0.784
[1] 0.7842
[1] 0.7844
[1] 0.7846
[1] 0.7848
[1] 0.785
[1] 0.7852
[1] 0.7854
[1] 0.7856
[1] 0.7858
[1] 0.786
[1] 0.7862
[1] 0.7864
[1] 0.7866
[1] 0.7868
[1] 0.787
[1] 0.7872
[1] 0.7874
[1] 0.7876
[1] 0.7878
[1] 0.788
[1] 0.7882
[1] 0.7884
[1] 0.7886
[1] 0.7888
[1] 0.789
[1] 0.7892
[1] 0.7894
[1] 0.7896
[1] 0.7898
[1] 0.79
[1] 0.7902
[1] 0.7904
[1] 0.7906
[1] 0.7908
[1] 0.791
[1] 0.7912
[1] 0.7914
[1] 0.7916
[1] 0.7918
[1] 0.792
[1] 0.7922
[1] 0.7924
[1] 0.7926
[1] 0.7928
[1] 0.793
[1] 0.7932
[1] 0.7934
[1] 0.7936
[1] 0.7938
[1] 0.794
[1] 0.7942
[1] 0.7944
[1] 0.7946
[1] 0.7948
[1] 0.795
[1] 0.7952
[1] 0.7954
[1] 0.7956
[1] 0.7958
[1] 0.796
[1] 0.7962
[1] 0.7964
[1] 0.7966
[1] 0.7968
[1] 0.797
[1] 0.7972
[1] 0.7974
[1] 0.7976
[1] 0.7978
[1] 0.798
[1] 0.7982
[1] 0.7984
[1] 0.7986
[1] 0.7988
[1] 0.799
[1] 0.7992
[1] 0.7994
[1] 0.7996
[1] 0.7998
[1] 0.8
[1] 0.8002
[1] 0.8004
[1] 0.8006
[1] 0.8008
[1] 0.801
[1] 0.8012
[1] 0.8014
[1] 0.8016
[1] 0.8018
[1] 0.802
[1] 0.8022
[1] 0.8024
[1] 0.8026
[1] 0.8028
[1] 0.803
[1] 0.8032
[1] 0.8034
[1] 0.8036
[1] 0.8038
[1] 0.804
[1] 0.8042
[1] 0.8044
[1] 0.8046
[1] 0.8048
[1] 0.805
[1] 0.8052
[1] 0.8054
[1] 0.8056
[1] 0.8058
[1] 0.806
[1] 0.8062
[1] 0.8064
[1] 0.8066
[1] 0.8068
[1] 0.807
[1] 0.8072
[1] 0.8074
[1] 0.8076
[1] 0.8078
[1] 0.808
[1] 0.8082
[1] 0.8084
[1] 0.8086
[1] 0.8088
[1] 0.809
[1] 0.8092
[1] 0.8094
[1] 0.8096
[1] 0.8098
[1] 0.81
[1] 0.8102
[1] 0.8104
[1] 0.8106
[1] 0.8108
[1] 0.811
[1] 0.8112
[1] 0.8114
[1] 0.8116
[1] 0.8118
[1] 0.812
[1] 0.8122
[1] 0.8124
[1] 0.8126
[1] 0.8128
[1] 0.813
[1] 0.8132
[1] 0.8134
[1] 0.8136
[1] 0.8138
[1] 0.814
[1] 0.8142
[1] 0.8144
[1] 0.8146
[1] 0.8148
[1] 0.815
[1] 0.8152
[1] 0.8154
[1] 0.8156
[1] 0.8158
[1] 0.816
[1] 0.8162
[1] 0.8164
[1] 0.8166
[1] 0.8168
[1] 0.817
[1] 0.8172
[1] 0.8174
[1] 0.8176
[1] 0.8178
[1] 0.818
[1] 0.8182
[1] 0.8184
[1] 0.8186
[1] 0.8188
[1] 0.819
[1] 0.8192
[1] 0.8194
[1] 0.8196
[1] 0.8198
[1] 0.82
[1] 0.8202
[1] 0.8204
[1] 0.8206
[1] 0.8208
[1] 0.821
[1] 0.8212
[1] 0.8214
[1] 0.8216
[1] 0.8218
[1] 0.822
[1] 0.8222
[1] 0.8224
[1] 0.8226
[1] 0.8228
[1] 0.823
[1] 0.8232
[1] 0.8234
[1] 0.8236
[1] 0.8238
[1] 0.824
[1] 0.8242
[1] 0.8244
[1] 0.8246
[1] 0.8248
[1] 0.825
[1] 0.8252
[1] 0.8254
[1] 0.8256
[1] 0.8258
[1] 0.826
[1] 0.8262
[1] 0.8264
[1] 0.8266
[1] 0.8268
[1] 0.827
[1] 0.8272
[1] 0.8274
[1] 0.8276
[1] 0.8278
[1] 0.828
[1] 0.8282
[1] 0.8284
[1] 0.8286
[1] 0.8288
[1] 0.829
[1] 0.8292
[1] 0.8294
[1] 0.8296
[1] 0.8298
[1] 0.83
[1] 0.8302
[1] 0.8304
[1] 0.8306
[1] 0.8308
[1] 0.831
[1] 0.8312
[1] 0.8314
[1] 0.8316
[1] 0.8318
[1] 0.832
[1] 0.8322
[1] 0.8324
[1] 0.8326
[1] 0.8328
[1] 0.833
[1] 0.8332
[1] 0.8334
[1] 0.8336
[1] 0.8338
[1] 0.834
[1] 0.8342
[1] 0.8344
[1] 0.8346
[1] 0.8348
[1] 0.835
[1] 0.8352
[1] 0.8354
[1] 0.8356
[1] 0.8358
[1] 0.836
[1] 0.8362
[1] 0.8364
[1] 0.8366
[1] 0.8368
[1] 0.837
[1] 0.8372
[1] 0.8374
[1] 0.8376
[1] 0.8378
[1] 0.838
[1] 0.8382
[1] 0.8384
[1] 0.8386
[1] 0.8388
[1] 0.839
[1] 0.8392
[1] 0.8394
[1] 0.8396
[1] 0.8398
[1] 0.84
[1] 0.8402
[1] 0.8404
[1] 0.8406
[1] 0.8408
[1] 0.841
[1] 0.8412
[1] 0.8414
[1] 0.8416
[1] 0.8418
[1] 0.842
[1] 0.8422
[1] 0.8424
[1] 0.8426
[1] 0.8428
[1] 0.843
[1] 0.8432
[1] 0.8434
[1] 0.8436
[1] 0.8438
[1] 0.844
[1] 0.8442
[1] 0.8444
[1] 0.8446
[1] 0.8448
[1] 0.845
[1] 0.8452
[1] 0.8454
[1] 0.8456
[1] 0.8458
[1] 0.846
[1] 0.8462
[1] 0.8464
[1] 0.8466
[1] 0.8468
[1] 0.847
[1] 0.8472
[1] 0.8474
[1] 0.8476
[1] 0.8478
[1] 0.848
[1] 0.8482
[1] 0.8484
[1] 0.8486
[1] 0.8488
[1] 0.849
[1] 0.8492
[1] 0.8494
[1] 0.8496
[1] 0.8498
[1] 0.85
[1] 0.8502
[1] 0.8504
[1] 0.8506
[1] 0.8508
[1] 0.851
[1] 0.8512
[1] 0.8514
[1] 0.8516
[1] 0.8518
[1] 0.852
[1] 0.8522
[1] 0.8524
[1] 0.8526
[1] 0.8528
[1] 0.853
[1] 0.8532
[1] 0.8534
[1] 0.8536
[1] 0.8538
[1] 0.854
[1] 0.8542
[1] 0.8544
[1] 0.8546
[1] 0.8548
[1] 0.855
[1] 0.8552
[1] 0.8554
[1] 0.8556
[1] 0.8558
[1] 0.856
[1] 0.8562
[1] 0.8564
[1] 0.8566
[1] 0.8568
[1] 0.857
[1] 0.8572
[1] 0.8574
[1] 0.8576
[1] 0.8578
[1] 0.858
[1] 0.8582
[1] 0.8584
[1] 0.8586
[1] 0.8588
[1] 0.859
[1] 0.8592
[1] 0.8594
[1] 0.8596
[1] 0.8598
[1] 0.86
[1] 0.8602
[1] 0.8604
[1] 0.8606
[1] 0.8608
[1] 0.861
[1] 0.8612
[1] 0.8614
[1] 0.8616
[1] 0.8618
[1] 0.862
[1] 0.8622
[1] 0.8624
[1] 0.8626
[1] 0.8628
[1] 0.863
[1] 0.8632
[1] 0.8634
[1] 0.8636
[1] 0.8638
[1] 0.864
[1] 0.8642
[1] 0.8644
[1] 0.8646
[1] 0.8648
[1] 0.865
[1] 0.8652
[1] 0.8654
[1] 0.8656
[1] 0.8658
[1] 0.866
[1] 0.8662
[1] 0.8664
[1] 0.8666
[1] 0.8668
[1] 0.867
[1] 0.8672
[1] 0.8674
[1] 0.8676
[1] 0.8678
[1] 0.868
[1] 0.8682
[1] 0.8684
[1] 0.8686
[1] 0.8688
[1] 0.869
[1] 0.8692
[1] 0.8694
[1] 0.8696
[1] 0.8698
[1] 0.87
[1] 0.8702
[1] 0.8704
[1] 0.8706
[1] 0.8708
[1] 0.871
[1] 0.8712
[1] 0.8714
[1] 0.8716
[1] 0.8718
[1] 0.872
[1] 0.8722
[1] 0.8724
[1] 0.8726
[1] 0.8728
[1] 0.873
[1] 0.8732
[1] 0.8734
[1] 0.8736
[1] 0.8738
[1] 0.874
[1] 0.8742
[1] 0.8744
[1] 0.8746
[1] 0.8748
[1] 0.875
[1] 0.8752
[1] 0.8754
[1] 0.8756
[1] 0.8758
[1] 0.876
[1] 0.8762
[1] 0.8764
[1] 0.8766
[1] 0.8768
[1] 0.877
[1] 0.8772
[1] 0.8774
[1] 0.8776
[1] 0.8778
[1] 0.878
[1] 0.8782
[1] 0.8784
[1] 0.8786
[1] 0.8788
[1] 0.879
[1] 0.8792
[1] 0.8794
[1] 0.8796
[1] 0.8798
[1] 0.88
[1] 0.8802
[1] 0.8804
[1] 0.8806
[1] 0.8808
[1] 0.881
[1] 0.8812
[1] 0.8814
[1] 0.8816
[1] 0.8818
[1] 0.882
[1] 0.8822
[1] 0.8824
[1] 0.8826
[1] 0.8828
[1] 0.883
[1] 0.8832
[1] 0.8834
[1] 0.8836
[1] 0.8838
[1] 0.884
[1] 0.8842
[1] 0.8844
[1] 0.8846
[1] 0.8848
[1] 0.885
[1] 0.8852
[1] 0.8854
[1] 0.8856
[1] 0.8858
[1] 0.886
[1] 0.8862
[1] 0.8864
[1] 0.8866
[1] 0.8868
[1] 0.887
[1] 0.8872
[1] 0.8874
[1] 0.8876
[1] 0.8878
[1] 0.888
[1] 0.8882
[1] 0.8884
[1] 0.8886
[1] 0.8888
[1] 0.889
[1] 0.8892
[1] 0.8894
[1] 0.8896
[1] 0.8898
[1] 0.89
[1] 0.8902
[1] 0.8904
[1] 0.8906
[1] 0.8908
[1] 0.891
[1] 0.8912
[1] 0.8914
[1] 0.8916
[1] 0.8918
[1] 0.892
[1] 0.8922
[1] 0.8924
[1] 0.8926
[1] 0.8928
[1] 0.893
[1] 0.8932
[1] 0.8934
[1] 0.8936
[1] 0.8938
[1] 0.894
[1] 0.8942
[1] 0.8944
[1] 0.8946
[1] 0.8948
[1] 0.895
[1] 0.8952
[1] 0.8954
[1] 0.8956
[1] 0.8958
[1] 0.896
[1] 0.8962
[1] 0.8964
[1] 0.8966
[1] 0.8968
[1] 0.897
[1] 0.8972
[1] 0.8974
[1] 0.8976
[1] 0.8978
[1] 0.898
[1] 0.8982
[1] 0.8984
[1] 0.8986
[1] 0.8988
[1] 0.899
[1] 0.8992
[1] 0.8994
[1] 0.8996
[1] 0.8998
[1] 0.9
[1] 0.9002
[1] 0.9004
[1] 0.9006
[1] 0.9008
[1] 0.901
[1] 0.9012
[1] 0.9014
[1] 0.9016
[1] 0.9018
[1] 0.902
[1] 0.9022
[1] 0.9024
[1] 0.9026
[1] 0.9028
[1] 0.903
[1] 0.9032
[1] 0.9034
[1] 0.9036
[1] 0.9038
[1] 0.904
[1] 0.9042
[1] 0.9044
[1] 0.9046
[1] 0.9048
[1] 0.905
[1] 0.9052
[1] 0.9054
[1] 0.9056
[1] 0.9058
[1] 0.906
[1] 0.9062
[1] 0.9064
[1] 0.9066
[1] 0.9068
[1] 0.907
[1] 0.9072
[1] 0.9074
[1] 0.9076
[1] 0.9078
[1] 0.908
[1] 0.9082
[1] 0.9084
[1] 0.9086
[1] 0.9088
[1] 0.909
[1] 0.9092
[1] 0.9094
[1] 0.9096
[1] 0.9098
[1] 0.91
[1] 0.9102
[1] 0.9104
[1] 0.9106
[1] 0.9108
[1] 0.911
[1] 0.9112
[1] 0.9114
[1] 0.9116
[1] 0.9118
[1] 0.912
[1] 0.9122
[1] 0.9124
[1] 0.9126
[1] 0.9128
[1] 0.913
[1] 0.9132
[1] 0.9134
[1] 0.9136
[1] 0.9138
[1] 0.914
[1] 0.9142
[1] 0.9144
[1] 0.9146
[1] 0.9148
[1] 0.915
[1] 0.9152
[1] 0.9154
[1] 0.9156
[1] 0.9158
[1] 0.916
[1] 0.9162
[1] 0.9164
[1] 0.9166
[1] 0.9168
[1] 0.917
[1] 0.9172
[1] 0.9174
[1] 0.9176
[1] 0.9178
[1] 0.918
[1] 0.9182
[1] 0.9184
[1] 0.9186
[1] 0.9188
[1] 0.919
[1] 0.9192
[1] 0.9194
[1] 0.9196
[1] 0.9198
[1] 0.92
[1] 0.9202
[1] 0.9204
[1] 0.9206
[1] 0.9208
[1] 0.921
[1] 0.9212
[1] 0.9214
[1] 0.9216
[1] 0.9218
[1] 0.922
[1] 0.9222
[1] 0.9224
[1] 0.9226
[1] 0.9228
[1] 0.923
[1] 0.9232
[1] 0.9234
[1] 0.9236
[1] 0.9238
[1] 0.924
[1] 0.9242
[1] 0.9244
[1] 0.9246
[1] 0.9248
[1] 0.925
[1] 0.9252
[1] 0.9254
[1] 0.9256
[1] 0.9258
[1] 0.926
[1] 0.9262
[1] 0.9264
[1] 0.9266
[1] 0.9268
[1] 0.927
[1] 0.9272
[1] 0.9274
[1] 0.9276
[1] 0.9278
[1] 0.928
[1] 0.9282
[1] 0.9284
[1] 0.9286
[1] 0.9288
[1] 0.929
[1] 0.9292
[1] 0.9294
[1] 0.9296
[1] 0.9298
[1] 0.93
[1] 0.9302
[1] 0.9304
[1] 0.9306
[1] 0.9308
[1] 0.931
[1] 0.9312
[1] 0.9314
[1] 0.9316
[1] 0.9318
[1] 0.932
[1] 0.9322
[1] 0.9324
[1] 0.9326
[1] 0.9328
[1] 0.933
[1] 0.9332
[1] 0.9334
[1] 0.9336
[1] 0.9338
[1] 0.934
[1] 0.9342
[1] 0.9344
[1] 0.9346
[1] 0.9348
[1] 0.935
[1] 0.9352
[1] 0.9354
[1] 0.9356
[1] 0.9358
[1] 0.936
[1] 0.9362
[1] 0.9364
[1] 0.9366
[1] 0.9368
[1] 0.937
[1] 0.9372
[1] 0.9374
[1] 0.9376
[1] 0.9378
[1] 0.938
[1] 0.9382
[1] 0.9384
[1] 0.9386
[1] 0.9388
[1] 0.939
[1] 0.9392
[1] 0.9394
[1] 0.9396
[1] 0.9398
[1] 0.94
[1] 0.9402
[1] 0.9404
[1] 0.9406
[1] 0.9408
[1] 0.941
[1] 0.9412
[1] 0.9414
[1] 0.9416
[1] 0.9418
[1] 0.942
[1] 0.9422
[1] 0.9424
[1] 0.9426
[1] 0.9428
[1] 0.943
[1] 0.9432
[1] 0.9434
[1] 0.9436
[1] 0.9438
[1] 0.944
[1] 0.9442
[1] 0.9444
[1] 0.9446
[1] 0.9448
[1] 0.945
[1] 0.9452
[1] 0.9454
[1] 0.9456
[1] 0.9458
[1] 0.946
[1] 0.9462
[1] 0.9464
[1] 0.9466
[1] 0.9468
[1] 0.947
[1] 0.9472
[1] 0.9474
[1] 0.9476
[1] 0.9478
[1] 0.948
[1] 0.9482
[1] 0.9484
[1] 0.9486
[1] 0.9488
[1] 0.949
[1] 0.9492
[1] 0.9494
[1] 0.9496
[1] 0.9498
[1] 0.95
[1] 0.9502
[1] 0.9504
[1] 0.9506
[1] 0.9508
[1] 0.951
[1] 0.9512
[1] 0.9514
[1] 0.9516
[1] 0.9518
[1] 0.952
[1] 0.9522
[1] 0.9524
[1] 0.9526
[1] 0.9528
[1] 0.953
[1] 0.9532
[1] 0.9534
[1] 0.9536
[1] 0.9538
[1] 0.954
[1] 0.9542
[1] 0.9544
[1] 0.9546
[1] 0.9548
[1] 0.955
[1] 0.9552
[1] 0.9554
[1] 0.9556
[1] 0.9558
[1] 0.956
[1] 0.9562
[1] 0.9564
[1] 0.9566
[1] 0.9568
[1] 0.957
[1] 0.9572
[1] 0.9574
[1] 0.9576
[1] 0.9578
[1] 0.958
[1] 0.9582
[1] 0.9584
[1] 0.9586
[1] 0.9588
[1] 0.959
[1] 0.9592
[1] 0.9594
[1] 0.9596
[1] 0.9598
[1] 0.96
[1] 0.9602
[1] 0.9604
[1] 0.9606
[1] 0.9608
[1] 0.961
[1] 0.9612
[1] 0.9614
[1] 0.9616
[1] 0.9618
[1] 0.962
[1] 0.9622
[1] 0.9624
[1] 0.9626
[1] 0.9628
[1] 0.963
[1] 0.9632
[1] 0.9634
[1] 0.9636
[1] 0.9638
[1] 0.964
[1] 0.9642
[1] 0.9644
[1] 0.9646
[1] 0.9648
[1] 0.965
[1] 0.9652
[1] 0.9654
[1] 0.9656
[1] 0.9658
[1] 0.966
[1] 0.9662
[1] 0.9664
[1] 0.9666
[1] 0.9668
[1] 0.967
[1] 0.9672
[1] 0.9674
[1] 0.9676
[1] 0.9678
[1] 0.968
[1] 0.9682
[1] 0.9684
[1] 0.9686
[1] 0.9688
[1] 0.969
[1] 0.9692
[1] 0.9694
[1] 0.9696
[1] 0.9698
[1] 0.97
[1] 0.9702
[1] 0.9704
[1] 0.9706
[1] 0.9708
[1] 0.971
[1] 0.9712
[1] 0.9714
[1] 0.9716
[1] 0.9718
[1] 0.972
[1] 0.9722
[1] 0.9724
[1] 0.9726
[1] 0.9728
[1] 0.973
[1] 0.9732
[1] 0.9734
[1] 0.9736
[1] 0.9738
[1] 0.974
[1] 0.9742
[1] 0.9744
[1] 0.9746
[1] 0.9748
[1] 0.975
[1] 0.9752
[1] 0.9754
[1] 0.9756
[1] 0.9758
[1] 0.976
[1] 0.9762
[1] 0.9764
[1] 0.9766
[1] 0.9768
[1] 0.977
[1] 0.9772
[1] 0.9774
[1] 0.9776
[1] 0.9778
[1] 0.978
[1] 0.9782
[1] 0.9784
[1] 0.9786
[1] 0.9788
[1] 0.979
[1] 0.9792
[1] 0.9794
[1] 0.9796
[1] 0.9798
[1] 0.98
[1] 0.9802
[1] 0.9804
[1] 0.9806
[1] 0.9808
[1] 0.981
[1] 0.9812
[1] 0.9814
[1] 0.9816
[1] 0.9818
[1] 0.982
[1] 0.9822
[1] 0.9824
[1] 0.9826
[1] 0.9828
[1] 0.983
[1] 0.9832
[1] 0.9834
[1] 0.9836
[1] 0.9838
[1] 0.984
[1] 0.9842
[1] 0.9844
[1] 0.9846
[1] 0.9848
[1] 0.985
[1] 0.9852
[1] 0.9854
[1] 0.9856
[1] 0.9858
[1] 0.986
[1] 0.9862
[1] 0.9864
[1] 0.9866
[1] 0.9868
[1] 0.987
[1] 0.9872
[1] 0.9874
[1] 0.9876
[1] 0.9878
[1] 0.988
[1] 0.9882
[1] 0.9884
[1] 0.9886
[1] 0.9888
[1] 0.989
[1] 0.9892
[1] 0.9894
[1] 0.9896
[1] 0.9898
[1] 0.99
[1] 0.9902
[1] 0.9904
[1] 0.9906
[1] 0.9908
[1] 0.991
[1] 0.9912
[1] 0.9914
[1] 0.9916
[1] 0.9918
[1] 0.992
[1] 0.9922
[1] 0.9924
[1] 0.9926
[1] 0.9928
[1] 0.993
[1] 0.9932
[1] 0.9934
[1] 0.9936
[1] 0.9938
[1] 0.994
[1] 0.9942
[1] 0.9944
[1] 0.9946
[1] 0.9948
[1] 0.995
[1] 0.9952
[1] 0.9954
[1] 0.9956
[1] 0.9958
[1] 0.996
[1] 0.9962
[1] 0.9964
[1] 0.9966
[1] 0.9968
[1] 0.997
[1] 0.9972
[1] 0.9974
[1] 0.9976
[1] 0.9978
[1] 0.998
[1] 0.9982
[1] 0.9984
[1] 0.9986
[1] 0.9988
[1] 0.999
[1] 0.9992
[1] 0.9994
[1] 0.9996
[1] 0.9998
[1] 1
if (sum(abs(Y) != abs(Re(Y))) == 0) {
  Y <- Re(Y)
}
nrow(input_var)
[1] 5000
nrow(Y)
[1] 5000
Y <- as.data.frame(Y); colnames(Y) <- c("eigenvalue2", "trace")
(data2 <- cbind(input_var, Y)) %>% head()
data2$p <- factor(data2$p)
ggplot(data2[30:nrow(data2),], mapping = aes(col = p, x = beta, y = eigenvalue2)) +
  geom_point(size = 0.5) +
  stat_smooth() +
  scale_size_continuous(range = c(0.1,1)) +
  labs(x = "Beta", y = TeX("Second Eigenvalue"), title = TeX("Second Eigenvalue given Matrix Density") ) +
  guides(col = guide_legend("Link Density"))  +
  theme_bw()

These look logarithmic, let’s do a log transform:


ggplot(data2[30:nrow(data2),], mapping = aes(col = p, x = log(beta), y = eigenvalue2)) +
  geom_point(size = 0.5) +
  stat_smooth() +
  scale_size_continuous(range = c(0.1,1)) +
  labs(x = "Beta", y = TeX("Second Eigenvalue"), title = TeX("Second Eigenvalue given Matrix Density") ) +
  guides(col = guide_legend("Link Density"))  +
  theme_bw()

Link density is a function of size so E2ln(Beta) is still pretty useful.

it seems that this is very nearly linear for low densities, the web would have a low density, for example ba graphs appear bounded above by 1/2:

n_vec <- 0:20
m <- 3
y <- c()
for (n in n_vec) {
    x <- sample_pa(n = n, power = 3, m = m) %>% 
      get.adjacency() %>% 
      mean()
    y <- c(y, x)
}
plot(n_vec, y)

length(n_vec)
[1] 21
length(y)
[1] 21

So considering small Densities is wise, let’s have a look.

So at this stage I would like to fit a model over what we have in the first plotly diagram, but, I’m going to look at ba graphs in case there is something more insightful.


ggplot(data2[30:nrow(data2),], mapping = aes(col = factor(p), x = log(beta*p), y = eigenvalue2)) +
  geom_point(size = 0.5) +
  stat_smooth() +
  scale_size_continuous(range = c(0.1,1)) +
  labs(x = "Beta", y = TeX("Second Eigenvalue"), title = TeX("Second Eigenvalue given Matrix Density") ) +
  
  guides(col = guide_legend("Link Density"))  +
  theme_bw()

It seems for very low densities that it is linear and then logarithmic for moderate densities. BA models have density~=m/n so for the BA model we’ll just use beta.

p       <- seq(from = 0.1, to = 0.95, length.out = 40)
size    <- seq(from = 100, to = 1000, length.out = 5)
beta <- seq(from = 1, to = 10, length.out = 1000)
beta <- 5
input_var <- expand.grid("p" = p, "beta" = beta, "size" = size)
input_var

nc <- length(random_graph(1, 1, 1))
Y <- matrix(ncol = nc, nrow = nrow(input_var))
for (i in 1:nrow(input_var)) {
  X <- as.vector(input_var[i,])
  Y[i,] <-  random_graph(X$p, X$beta, X$size)
  print(i/nrow(input_var))
}
if (sum(abs(Y) != abs(Re(Y))) == 0) {
  Y <- Re(Y)
}
nrow(input_var)
nrow(Y)
Y <- as.data.frame(Y); colnames(Y) <- c("eigenvalue2", "trace")
(data2 <- cbind(input_var, Y)) %>% head()
data2$p <- factor(data2$p)



ggplot(data2, mapping = aes(col = factor(p), x = log(size), y = eigenvalue2)) +
  geom_point(size = 0.5) +
#  stat_smooth() +
#  scale_size_continuous(range = c(0.1,1)) +
  labs(x = "P", y = TeX("Second Eigenvalue"), title = TeX("Second Eigenvalue given Matrix Density") ) +
  guides(col = guide_legend("Size"))  +
  theme_bw()
LS0tCnRpdGxlOiAiSW52ZXN0aWdhdGluZyBTZWNvbmQgRWlnZW52YWx1ZSIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKIyBJbnZlc3RpZ2F0aW5nIFNlY29uZCBFaWdlbnZhbHVlCgpgYGB7cn0KICBpZiAocmVxdWlyZSgiUGFnZVJhbmsiKSkgewogICAgICBsaWJyYXJ5KFBhZ2VSYW5rKQogICAgfWVsc2V7CiAgICAgIGRldnRvb2xzOjppbnN0YWxsX2dpdGh1YigicnlhbmdyZWVudXAvUGFnZVJhbmsiKQogICAgICBsaWJyYXJ5KFBhZ2VSYW5rKQogICAgfQoKICBsaWJyYXJ5KHBhY21hbikKICBwYWNtYW46OnBfbG9hZChQYWdlUmFuaywgZGV2dG9vbHMsIE1hdHJpeCwgaWdyYXBoLCBtaXNlLCB0aWR5dmVyc2UsIHJnbCwgbGF0ZXgyZXhwKQojICBtaXNlKCkKYGBgCgoKCiMjIExvb2tpbmcgYXQgRGVuc2l0eQoKIyMjIENvbnN0YW50cwoKRGVmaW5lIHNvbWUgY29uc3RhbnRzCgpgYGB7cn0KcCAgICA8LSBzZXEoZnJvbSA9IDAuMDEsIHRvID0gMC45OSwgbGVuZ3RoLm91dCA9IDEwKQpiZXRhIDwtIHNlcShmcm9tID0gMSwgdG8gPSAyMCwgbGVuZ3RoLm91dCA9IDQpCnN6IDwtIHNlcShmcm9tID0gMTAwLCB0byA9IDEwLCBsZW5ndGgub3V0ID0gNSkKaW5wdXRfdmFyIDwtIGV4cGFuZC5ncmlkKCJwIiA9IHAsICJiZXRhIiA9IGJldGEsICJzaXplIiA9IHN6KQppbnB1dF92YXIKYGBgCgojIyMgRnVuY3Rpb24gdG8gQnVpbGQgR3JhcGgKCmBgYHtyfQpyYW5kb21fZ3JhcGggPC0gZnVuY3Rpb24ocCwgYmV0YSwgc2l6ZSkgewogICAgICBnMSA8LSBpZ3JhcGg6OmVyZG9zLnJlbnlpLmdhbWUobiA9IHNpemUsIHApCiAgICAgIEEgPC0gaWdyYXBoOjpnZXQuYWRqYWNlbmN5KGcxKSAjIFJvdyB0byBjb2x1bW4KICAgICAgQSA8LSBNYXRyaXg6OnQoQSkKCiAgICAgIEFfZGVucyA8LSBtZWFuKEEpCiAgICAgIFQgICAgICA8LSBQYWdlUmFuazo6cG93ZXJfd2Fsa19wcm9iX3RyYW5zKEEsIGJldGEgPSBiZXRhKQogICAgICB0ciAgICAgPC0gc3VtKGRpYWcoVCkpCiAgICAgIGUyICAgICA8LSBlaWdlbihULCBvbmx5LnZhbHVlcyA9IFRSVUUpJHZhbHVlc1syXSAjIFIgb3JkZXJzIGJ5IGRlc2NlbmRpbmcgbWFnbml0dWRlCiAgICAgIHJldHVybihjKGFicyhlMiksIG1lYW4oQSksIHRyKSkKfQpgYGAKCiMjIyBSZXR1cm4gcmVzdWx0cwoKTWFwIHRoZSBmdW5jdGlvbgoKYGBge3J9Cm5jIDwtIGxlbmd0aChyYW5kb21fZ3JhcGgoMSwgMSwgMSkpClkgPC0gbWF0cml4KG5jb2wgPSBuYywgbnJvdyA9IG5yb3coaW5wdXRfdmFyKSkKZm9yIChpIGluIDE6bnJvdyhpbnB1dF92YXIpKSB7CiAgWCA8LSBhcy52ZWN0b3IoaW5wdXRfdmFyW2ksXSkKICBZW2ksXSA8LSAgcmFuZG9tX2dyYXBoKFgkcCwgWCRiZXRhLCBYJHNpemUpCiAgcHJpbnQoaS9ucm93KGlucHV0X3ZhcikpCn0KaWYgKHN1bShhYnMoWSkgIT0gYWJzKFJlKFkpKSkgPT0gMCkgewogIFkgPC0gUmUoWSkKfQpucm93KGlucHV0X3ZhcikKbnJvdyhZKQpZIDwtIGFzLmRhdGEuZnJhbWUoWSk7IGNvbG5hbWVzKFkpIDwtIGMoImVpZ2VudmFsdWUyIiwgIkFfZGVucyIsICJ0cmFjZSIpCihkYXRhIDwtIGNiaW5kKGlucHV0X3ZhciwgWSkpICU+JSBoZWFkKCkKYGBgCgoKIyMjIFBsb3QgUmVzdWx0cwoKCmBgYHtyfQpwYWlycyhkYXRhKQpjb3IoZGF0YSkKbGlicmFyeShjb3JycGxvdCkKY29ybWF0ID0gY29yKGRhdGEsIG1ldGhvZCA9ICdzcGVhcm1hbicpCmNvcnJwbG90KGNvcm1hdCwgbWV0aG9kID0gImVsbGlwc2UiLCB0eXBlID0gImxvd2VyIikKbmFtZXMoZGF0YSkKYGBgCkxldCdzIGxvb2sgYXQgYSAzZCBPdXRwdXQKCmBgYHtyfQojIHBsb3QzZChkYXRhJGJldGEsIGRhdGEkQV9kZW5zLCBkYXRhJGVpZ2VudmFsdWUyLCBzaXplID0gMikKYGBgCgoKCmBgYHtyfQoKbmFtZXMoZGF0YSkKbGlicmFyeShwbG90bHkpCgpkIDwtIGRhdGEKIyBkJGJldGEgPC0gbG9nKGQkYmV0YSkKCmZpZyA8LSBwbG90X2x5KGQsIHggPSB+QV9kZW5zLCB5ID0gfmJldGEsIHogPSB+ZWlnZW52YWx1ZTIpCmZpZyA8LSBmaWcgJT4lIGFkZF9tYXJrZXJzKHNpemUgPSAxKQpmaWcgPC0gZmlnICU+JSBsYXlvdXQoc2NlbmUgPSBsaXN0KHhheGlzID0gbGlzdCh0aXRsZSA9ICdEZW5zaXR5JyksCiAgICAgICAgICAgICAgICAgICAgIHlheGlzID0gbGlzdCh0aXRsZSA9ICdCZXRhJyksCiAgICAgICAgICAgICAgICAgICAgIHpheGlzID0gbGlzdCh0aXRsZSA9ICdFMicpKSkKCmZpZwoKYGBgCgpgYGB7cn0KCm5hbWVzKGRhdGEpCmxpYnJhcnkocGxvdGx5KQoKI2QgPC0gZGF0YVtzYW1wbGUoMTpucm93KGRhdGEpLCAxMDAwKSxdCmQgPC0gZGF0YQoKZmlnIDwtIHBsb3RfbHkoZCwgeCA9IH5BX2RlbnMsIHkgPSB+dHJhY2UsIHogPSB+ZWlnZW52YWx1ZTIpCmZpZyA8LSBmaWcgJT4lIGFkZF9tYXJrZXJzKHNpemUgPSAxKQpmaWcgPC0gZmlnICU+JSBsYXlvdXQoc2NlbmUgPSBsaXN0KHhheGlzID0gbGlzdCh0aXRsZSA9ICdEZW5zaXR5JyksCiAgICAgICAgICAgICAgICAgICAgIHlheGlzID0gbGlzdCh0aXRsZSA9ICd0cmFjZScpLAogICAgICAgICAgICAgICAgICAgICB6YXhpcyA9IGxpc3QodGl0bGUgPSAnRTInKSkpCgpmaWcKCmBgYAoKCkNsZWFybHkgSSBzaG91bGQgYmUgYWJsZSB0byBtb2RlbCB0aGlzLgoKTGV0J3MgbG9vayBhdCBkZW5zaXR5LCBpdCdzIGNvbnRpbnVvdXMgc28gSSBzaG91bGQgYmUgYWJsZSB0byAgdGFrZSBzcGxpY2VzIG9mIEFfZGVucy4KCkFfZGVucyBpcyBkZXBlbmRlbnQgb24gcCwgc28gSSdsbCBqdXN0IHVzZSBkaWZmZXJlbnQgcHMKCgpgYGB7cn0KcCAgICA8LSBzZXEoZnJvbSA9IDAuMSwgdG8gPSAwLjk1LCBsZW5ndGgub3V0ID0gNSkKYmV0YSA8LSBzZXEoZnJvbSA9IDEsIHRvID0gMTAsIGxlbmd0aC5vdXQgPSAxMDAwKQpzaXplID0gMTAwCmlucHV0X3ZhciA8LSBleHBhbmQuZ3JpZCgicCIgPSBwLCAiYmV0YSIgPSBiZXRhLCAic2l6ZSIgPSBzaXplKQppbnB1dF92YXIKCm5jIDwtIGxlbmd0aChyYW5kb21fZ3JhcGgoMSwgMSwgMSkpClkgPC0gbWF0cml4KG5jb2wgPSBuYywgbnJvdyA9IG5yb3coaW5wdXRfdmFyKSkKZm9yIChpIGluIDE6bnJvdyhpbnB1dF92YXIpKSB7CiAgWCA8LSBhcy52ZWN0b3IoaW5wdXRfdmFyW2ksXSkKICBZW2ksXSA8LSAgcmFuZG9tX2dyYXBoKFgkcCwgWCRiZXRhLCBYJHNpemUpCiAgcHJpbnQoaS9ucm93KGlucHV0X3ZhcikpCn0KaWYgKHN1bShhYnMoWSkgIT0gYWJzKFJlKFkpKSkgPT0gMCkgewogIFkgPC0gUmUoWSkKfQpucm93KGlucHV0X3ZhcikKbnJvdyhZKQpZIDwtIGFzLmRhdGEuZnJhbWUoWSk7IGNvbG5hbWVzKFkpIDwtIGMoImVpZ2VudmFsdWUyIiwgInRyYWNlIikKKGRhdGEyIDwtIGNiaW5kKGlucHV0X3ZhciwgWSkpICU+JSBoZWFkKCkKZGF0YTIkcCA8LSBmYWN0b3IoZGF0YTIkcCkKCmBgYAoKCmBgYHtyfQpnZ3Bsb3QoZGF0YTJbMzA6bnJvdyhkYXRhMiksXSwgbWFwcGluZyA9IGFlcyhjb2wgPSBwLCB4ID0gYmV0YSwgeSA9IGVpZ2VudmFsdWUyKSkgKwogIGdlb21fcG9pbnQoc2l6ZSA9IDAuNSkgKwogIHN0YXRfc21vb3RoKCkgKwogIHNjYWxlX3NpemVfY29udGludW91cyhyYW5nZSA9IGMoMC4xLDEpKSArCiAgbGFicyh4ID0gIkJldGEiLCB5ID0gVGVYKCJTZWNvbmQgRWlnZW52YWx1ZSIpLCB0aXRsZSA9IFRlWCgiU2Vjb25kIEVpZ2VudmFsdWUgZ2l2ZW4gTWF0cml4IERlbnNpdHkiKSApICsKICBndWlkZXMoY29sID0gZ3VpZGVfbGVnZW5kKCJMaW5rIERlbnNpdHkiKSkgICsKICB0aGVtZV9idygpCmBgYAoKVGhlc2UgbG9vayBsb2dhcml0aG1pYywgbGV0J3MgZG8gYSBsb2cgdHJhbnNmb3JtOgoKYGBge3J9CgpnZ3Bsb3QoZGF0YTJbMzA6bnJvdyhkYXRhMiksXSwgbWFwcGluZyA9IGFlcyhjb2wgPSBwLCB4ID0gbG9nKGJldGEpLCB5ID0gZWlnZW52YWx1ZTIpKSArCiAgZ2VvbV9wb2ludChzaXplID0gMC41KSArCiAgc3RhdF9zbW9vdGgoKSArCiAgc2NhbGVfc2l6ZV9jb250aW51b3VzKHJhbmdlID0gYygwLjEsMSkpICsKICBsYWJzKHggPSAiQmV0YSIsIHkgPSBUZVgoIlNlY29uZCBFaWdlbnZhbHVlIiksIHRpdGxlID0gVGVYKCJTZWNvbmQgRWlnZW52YWx1ZSBnaXZlbiBNYXRyaXggRGVuc2l0eSIpICkgKwogIGd1aWRlcyhjb2wgPSBndWlkZV9sZWdlbmQoIkxpbmsgRGVuc2l0eSIpKSAgKwogIHRoZW1lX2J3KCkKYGBgCkxpbmsgZGVuc2l0eSBpcyBhIGZ1bmN0aW9uIG9mIHNpemUgc28gRTJccHJvcHRvIGxuKEJldGEpIGlzIHN0aWxsIHByZXR0eSB1c2VmdWwuCgppdCBzZWVtcyB0aGF0IHRoaXMgaXMgdmVyeSBuZWFybHkgbGluZWFyIGZvciBsb3cgZGVuc2l0aWVzLCB0aGUgd2ViIHdvdWxkIGhhdmUgYSBsb3cgZGVuc2l0eSwgZm9yIGV4YW1wbGUgYmEgZ3JhcGhzIGFwcGVhciBib3VuZGVkIGFib3ZlIGJ5IDEvMjoKCmBgYHtyfQpuX3ZlYyA8LSAwOjIwCm0gPC0gMwp5IDwtIGMoKQpmb3IgKG4gaW4gbl92ZWMpIHsKICAgIHggPC0gc2FtcGxlX3BhKG4gPSBuLCBwb3dlciA9IDMsIG0gPSBtKSAlPiUgCiAgICAgIGdldC5hZGphY2VuY3koKSAlPiUgCiAgICAgIG1lYW4oKQogICAgeSA8LSBjKHksIHgpCn0KcGxvdChuX3ZlYywgeSkKbGVuZ3RoKG5fdmVjKQpsZW5ndGgoeSkKYGBgCgpTbyBjb25zaWRlcmluZyBzbWFsbCBEZW5zaXRpZXMgaXMgd2lzZSwgbGV0J3MgaGF2ZSBhIGxvb2suCgpTbyBhdCB0aGlzIHN0YWdlIEkgd291bGQgbGlrZSB0byBmaXQgYSBtb2RlbCBvdmVyIHdoYXQgd2UgaGF2ZSBpbiB0aGUgZmlyc3QgcGxvdGx5IGRpYWdyYW0sIGJ1dCwgSSdtIGdvaW5nIHRvIGxvb2sgYXQgYmEgZ3JhcGhzIGluIGNhc2UgdGhlcmUgaXMgc29tZXRoaW5nIG1vcmUgaW5zaWdodGZ1bC4KCmBgYHtyfQoKZ2dwbG90KGRhdGEyWzMwOm5yb3coZGF0YTIpLF0sIG1hcHBpbmcgPSBhZXMoY29sID0gZmFjdG9yKHApLCB4ID0gbG9nKGJldGEqcCksIHkgPSBlaWdlbnZhbHVlMikpICsKICBnZW9tX3BvaW50KHNpemUgPSAwLjUpICsKICBzdGF0X3Ntb290aCgpICsKICBzY2FsZV9zaXplX2NvbnRpbnVvdXMocmFuZ2UgPSBjKDAuMSwxKSkgKwogIGxhYnMoeCA9ICJCZXRhIiwgeSA9IFRlWCgiU2Vjb25kIEVpZ2VudmFsdWUiKSwgdGl0bGUgPSBUZVgoIlNlY29uZCBFaWdlbnZhbHVlIGdpdmVuIE1hdHJpeCBEZW5zaXR5IikgKSArCiAgCiAgZ3VpZGVzKGNvbCA9IGd1aWRlX2xlZ2VuZCgiTGluayBEZW5zaXR5IikpICArCiAgdGhlbWVfYncoKQpgYGAKCgoKSXQgc2VlbXMgZm9yIHZlcnkgbG93IGRlbnNpdGllcyB0aGF0IGl0IGlzIGxpbmVhciBhbmQgdGhlbiBsb2dhcml0aG1pYyBmb3IKbW9kZXJhdGUgZGVuc2l0aWVzLiBCQSBtb2RlbHMgaGF2ZSBkZW5zaXR5fj1tL24gc28gZm9yIHRoZSBCQSBtb2RlbCB3ZSdsbApqdXN0IHVzZSBiZXRhLgoKCmBgYHtyfQoKcCAgICAgICA8LSBzZXEoZnJvbSA9IDAuMSwgdG8gPSAwLjk1LCBsZW5ndGgub3V0ID0gNSkKc2l6ZSAgICA8LSBzZXEoZnJvbSA9IDEwMCwgdG8gPSAxMDAwLCBsZW5ndGgub3V0ID0gMzApCmJldGEgPC0gc2VxKGZyb20gPSAxLCB0byA9IDEwLCBsZW5ndGgub3V0ID0gMTAwMCkKYmV0YSA8LSA1CmlucHV0X3ZhciA8LSBleHBhbmQuZ3JpZCgicCIgPSBwLCAiYmV0YSIgPSBiZXRhLCAic2l6ZSIgPSBzaXplKQppbnB1dF92YXIKCm5jIDwtIGxlbmd0aChyYW5kb21fZ3JhcGgoMSwgMSwgMSkpClkgPC0gbWF0cml4KG5jb2wgPSBuYywgbnJvdyA9IG5yb3coaW5wdXRfdmFyKSkKZm9yIChpIGluIDE6bnJvdyhpbnB1dF92YXIpKSB7CiAgWCA8LSBhcy52ZWN0b3IoaW5wdXRfdmFyW2ksXSkKICBZW2ksXSA8LSAgcmFuZG9tX2dyYXBoKFgkcCwgWCRiZXRhLCBYJHNpemUpCiAgcHJpbnQoaS9ucm93KGlucHV0X3ZhcikpCn0KaWYgKHN1bShhYnMoWSkgIT0gYWJzKFJlKFkpKSkgPT0gMCkgewogIFkgPC0gUmUoWSkKfQpucm93KGlucHV0X3ZhcikKbnJvdyhZKQpZIDwtIGFzLmRhdGEuZnJhbWUoWSk7IGNvbG5hbWVzKFkpIDwtIGMoImVpZ2VudmFsdWUyIiwgInRyYWNlIikKKGRhdGEyIDwtIGNiaW5kKGlucHV0X3ZhciwgWSkpICU+JSBoZWFkKCkKZGF0YTIkcCA8LSBmYWN0b3IoZGF0YTIkcCkKCgoKZ2dwbG90KGRhdGEyLCBtYXBwaW5nID0gYWVzKGNvbCA9IGZhY3RvcihwKSwgeCA9IGxvZyhzaXplKSwgeSA9IGVpZ2VudmFsdWUyKSkgKwogIGdlb21fcG9pbnQoc2l6ZSA9IDAuNSkgKwojICBzdGF0X3Ntb290aCgpICsKIyAgc2NhbGVfc2l6ZV9jb250aW51b3VzKHJhbmdlID0gYygwLjEsMSkpICsKICBsYWJzKHggPSAic2l6ZSIsIHkgPSBUZVgoIlNlY29uZCBFaWdlbnZhbHVlIiksIHRpdGxlID0gVGVYKCJTZWNvbmQgRWlnZW52YWx1ZSBnaXZlbiBNYXRyaXggRGVuc2l0eSIpICkgKwogIGd1aWRlcyhjb2wgPSBndWlkZV9sZWdlbmQoIkxpbmsgRGVuc2l0eSIpKSAgKwogIHRoZW1lX2J3KCkKYGBgCgoKCgpgYGB7cn0KcCAgICAgICA8LSBzZXEoZnJvbSA9IDAuMSwgdG8gPSAwLjk1LCBsZW5ndGgub3V0ID0gNDApCnNpemUgICAgPC0gc2VxKGZyb20gPSAxMDAsIHRvID0gMTAwMCwgbGVuZ3RoLm91dCA9IDUpCmJldGEgPC0gc2VxKGZyb20gPSAxLCB0byA9IDEwLCBsZW5ndGgub3V0ID0gMTAwMCkKYmV0YSA8LSA1CmlucHV0X3ZhciA8LSBleHBhbmQuZ3JpZCgicCIgPSBwLCAiYmV0YSIgPSBiZXRhLCAic2l6ZSIgPSBzaXplKQppbnB1dF92YXIKCm5jIDwtIGxlbmd0aChyYW5kb21fZ3JhcGgoMSwgMSwgMSkpClkgPC0gbWF0cml4KG5jb2wgPSBuYywgbnJvdyA9IG5yb3coaW5wdXRfdmFyKSkKZm9yIChpIGluIDE6bnJvdyhpbnB1dF92YXIpKSB7CiAgWCA8LSBhcy52ZWN0b3IoaW5wdXRfdmFyW2ksXSkKICBZW2ksXSA8LSAgcmFuZG9tX2dyYXBoKFgkcCwgWCRiZXRhLCBYJHNpemUpCiAgcHJpbnQoaS9ucm93KGlucHV0X3ZhcikpCn0KaWYgKHN1bShhYnMoWSkgIT0gYWJzKFJlKFkpKSkgPT0gMCkgewogIFkgPC0gUmUoWSkKfQpucm93KGlucHV0X3ZhcikKbnJvdyhZKQpZIDwtIGFzLmRhdGEuZnJhbWUoWSk7IGNvbG5hbWVzKFkpIDwtIGMoImVpZ2VudmFsdWUyIiwgInRyYWNlIikKKGRhdGEyIDwtIGNiaW5kKGlucHV0X3ZhciwgWSkpICU+JSBoZWFkKCkKZGF0YTIkcCA8LSBmYWN0b3IoZGF0YTIkcCkKCgoKZ2dwbG90KGRhdGEyLCBtYXBwaW5nID0gYWVzKGNvbCA9IGZhY3RvcihwKSwgeCA9IGxvZyhzaXplKSwgeSA9IGVpZ2VudmFsdWUyKSkgKwogIGdlb21fcG9pbnQoc2l6ZSA9IDAuNSkgKwojICBzdGF0X3Ntb290aCgpICsKIyAgc2NhbGVfc2l6ZV9jb250aW51b3VzKHJhbmdlID0gYygwLjEsMSkpICsKICBsYWJzKHggPSAiUCIsIHkgPSBUZVgoIlNlY29uZCBFaWdlbnZhbHVlIiksIHRpdGxlID0gVGVYKCJTZWNvbmQgRWlnZW52YWx1ZSBnaXZlbiBNYXRyaXggRGVuc2l0eSIpICkgKwogIGd1aWRlcyhjb2wgPSBndWlkZV9sZWdlbmQoIlNpemUiKSkgICsKICB0aGVtZV9idygpCmBgYAoK